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ABSTRACT 


A  one-dimensional  mathematical  dynamic  simulation  model  for  predicting  the 
transient  behavior  of  air  induction  systems  for  advanced  tactical  aircraft  operating  in 
the  supersonic  portion  of  the  flight  spectrum  is  discussed.  As  a  part  of  the  simula¬ 
tion  model,  control  system  logic  is  included  to  provide  control  system  design  criteria 
necessary  to  maintain  air  induction  system  response  to  selected  input  disturbances 
within  prescribed  limits.  The  model  is  applicable  to  both  external  and  mixed  com¬ 
pression  inlets  operating  at  supersonic  conditions.  The  "lumped  parameter"  concept 
is  used  in  simulating  the  dynamic  response  of  the  subsonic  duct  downstream  of  the 
terminal  shock,  with  the  duct  divided  into  three  lumps.  CDC  6600  digital  computer 
programs  have  been  generated  for  each  of  three  sub-models  (i.e.,  Standard,  Buzz, 
and  Hammershock  models) .  Numerous  sample  simulation  runs  are  presented  and 
discussed  in  detail.  Limited  comparisons  of  simulation  results  to  available  test  data 
show  good  agreement. 
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SECTION  1 


INTRODUCTION 

The  basic  requirement  for  any  air  induction  system  is  to  provide  the  engine  with 
air  r.t  high  pressure  recovery  and  low  pressure  distortion  at  all  flight  conditions .  The 
increased  mission  requirements  for  future  highly  maneuverable  advanced  supersonic 
tactical  aircraft  have  made  the  task  of  designing  air  induction  systems  for  these  air¬ 
craft  extremely  complex.  In  addition  to  the  steady  state  requirements  for  high  effi¬ 
ciency  at  subsonic  and  supersonic  cruise  Mach  numbers,  the  air  induction  systems  for 
these  aircraft  must  be  designed  to  provide  stable  operation  of  the  propulsion  system 
during  transient  conditions  resulting  from  rapid  changes  in  engine  thrust,  flight  alti¬ 
tude,  and  aircraft  attitude  during  transonic  and  supersonic  maneuvering.  Thus,  pre¬ 
diction  of  the  transient  response  of  air  induction  systems  to  externally  or  internally 
generated  disturbances  whose  combined  rates  and  amplitudes  lead  to  non-equilibrium 
flow  phenomena  is  highly  desirable  during  preliminary  design  in  order  to  provide 
early  assessment  of  potential  problems  affecting  propulsion  system  stability. 

Reported  herein  is  a  one-dimensional  mathematical  dynamic  simulation  model 
for  predicting  the  transient  behavior  of  air  induction  systems  for  advanced  tactical 
aircraft  operating  in  the  supersonic  portion  of  the  flight  spectrum .  As  a  part  of  the 
simulation  model,  control  system  logic  is  included  to  provide  control  system  design 
criteria  necessary  to  maintain  air  induction  system  response  to  selected  input  dis¬ 
turbances  within  prescribed  limits.  The  model  is  applicable  to  both  external  and 
mixed  compression  inlets  operating  at  supersonic  conditions .  The  "lumped  parame¬ 
ter"  concept  is  used  in  simulating  the  dynamic  response  of  the  subsonic  duct  down¬ 
stream  of  the  terminal  shock,  with  the  duct  divided  into  three  lumps.  The  total 
dynamic  simulation  model  consists  of  three  sub-models  designated  as  follows: 

Model  1  -  Standard  Model 

Model  2  -  Buzz  Model 

Model  3  -  Hammer  shock  Model 
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The  decks  for  these  models  were  initially  developed  on  Northrop's  IBM  360  computer 
using  IBM  "Continuous  Systems  Modeling  Program"  (CSMP)  language  and  have  subse¬ 
quently  been  converted  to  a  similar  language  from  CDC  called  "Continuous  Systems 
Simulation  Language  (CSSL  3)"  for  use  on  the  Air  Force  CDC  6600  system.  A  user's 
manual  is  available  (Reference  1)  for  the  CDC  6600  (CSSL  3)  version  of  the  dynamic 
simulation  program.  It  describes  the  decks  in  detail  and  discusses  the  procedure  for 
operating  them . 

The  results  presented  herein  begin  with  a  brief  review  of  input  transients  likely 
to  be  encountered  by  advanced  tactical  aircraft.  Following  ihis  is  a  detailed  descrip¬ 
tion  of  the  theoretical  approach  used  in  developing  the  simulation  models,  with  sep¬ 
arate  derivations  of  key  equations  given  in  the  appendices.  Next,  numerous  sample 
simulation  runs  are  presented  and  discussed  in  detail .  Finally,  hammershock  simu¬ 
lation  results  are  compared  with  available  test  data . 
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SECTION  II 


AIK  INDUCTION  SYSTEM  THANS1E NTS 

Transient  disturbances  contributing  to  propulsion  system  instability  have  been 
well  documented  by  various  sources  (e.g.,  References  2  and  3).  These  disturbances 
may  be  generally  categorized  as  follows: 

1 .  External  Disturbances  -  Transient  disturbances  originating  external  to  the 
propulsion  system  such  as  wind  gusts,  clear  air  turbulence,  passing  air¬ 
craft  wakes,  aircraft  maneuvering,  hot  gas  ingestion  due  to  missile  and  gun 
firing,  etc. 

2.  Internal  Disturbances  -  Transient  disturbances  that  originate  within  the 
propulsion  system  such  as  sudden  changes  in  engine  airflow  demand  caused 
by  throttle  burst  or  chop,  afterburner  light-off  or  blowout,  engine  stall, 
etc. 

The  range  in  rates  and  amplitudes  of  input  disturbances  likely  to  be  encountered 
by  advanced  tactical  aircraft  is  wide.  Of  interest  with  respect  to  air  induction 
dynamic  simulation  are  those  transients  whose  combined  rates  and  amplitudes  are 
sufficiently  high  to  cause  significant  deviation  in  the  instantaneous  thermodynamic 
properties  within  the  air  induction  system  from  the  quasi-steady  state  properties 
associated  with  slower  transients.  For  example,  a  typical  engine  throttle  burst 
occurring  over  a  period  of  a  few  seconds,  while  seemingly  being  a  "transient"  situa¬ 
tion,  is  in  fact  essentially  a  quasi-steady  state  process  with  respect  to  the  response 
of  the  air  induction  system.  That  is,  at  each  instant  in  time  during  such  a  transient, 
the  airflow  within  the  induction  system  is  in  equilibrium  with  respect  to  the  instantane¬ 
ous  engine  demand. 

On  the  other  hand,  more  rapid  input  disturbances  (some  examples  of  which  are 
shown  in  Figure  1)  can  cause  significant  deviations  in  the  thermodynamic  properties 
within  the  air  induction  system  from  steady  state  conditions.  The  magnitude  of  devi¬ 
ation  from  steady  state  is  also  a  function  of  the  length/volume  characteristics  of  the 
air  induction  system.  For  these  cases,  non-equilibrium  effects  occur  due  to  the 
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Figure  1.  Typical  Air  Induction  System  Input  Transients 


inertial  and  capacitance  characteristics  of  air  within  the  induction  system .  The  result 
is  an  instantaneous  imbalance  of  airflows  and  pressures  within  the  duct  over  the  dura¬ 
tion  of  the  input  disturbance.  This  can  lead  to  propulsion  system  instability  due  to 
exposure  of  the  engine  compressor  face  to  rapid  rates  of  pressure  rise  and  decay 
and/or  adverse  instantaneous  pressure  distortion  and  turbulence  patterns .  It  can  also 
lead  to  excessive  pressure  levels  within  the  duct  (e.g.,  rapid  cut-off  of  engine  airflow 
can  result  in  instantaneous  duct  pressure  levels  in  excess  of  freestream  stagnation 
pressure) . 


SECTION  III 


SIMULATION  MODELS  -  THEORETICAL  APPROACH 

Dynamic  mathematical  models  are  currently  available  which  simulate  air  induc¬ 
tion  system  response  to  given  input  transients  using  various  calculation  approaches  in 
modeling  the  portion  of  the  subsonic  duct  downstream  of  the  normal  shock  (NS). 
Reference  4  uses  a  Helmholtz  resonator  spring-mass  concept  when  calculating  condi¬ 
tions  within  the  subsonic  duct.  Reference  5  uses  the  artificial  viscosity  technique  in 
calculating  ihe  effect  of  large  amplitude  transients  (i.e.,  compressor  surge).  Wasser- 
bauer  and  Willoh  (Reference  6)  use  a  linearized  small  perturbation  method  to  predict 
inlet  response,  while  Marvin  (Reference  2)  uses  the  'lumped  parameter'  concept  in 
simulating  air  induction  system  response  to  various  typical  transients. 

After  a  thorough  technical  review  of  the  various  available  models  and  an  assess¬ 
ment  of  their  suitability  for  achieving  the  objectives  of  Contract  F33615-69-C-1699, 
"Supersonic  Inlet  Investigation,  "  the  'lumped  parameter'  concept  developed  in  Refer¬ 
ence  2  was  selected,  based  on  the  merits  of  its  thermodynamic  formulation,  and  also 
because  the  range  of  models  developed  and  the  types  of  air  induction  systems  simu¬ 
lated  were  the  most  suitable  for  achieving  the  contract  objectives.  This  concept 
(which  is  explained  in  detail  in  this  section)  divides  the  subsonic  duct  into  a  number  of 
lumped  volumes.  Models  that  treat  the  whole  subsonic  duct  as  a  single  lumped  volume 
may  be  sufficiently  accurate  for  fairly  short  air  induction  systems  (podded  or  wing 
mounted  turbofan  installations),  but  their  accuracy  deteriorates  for  inlets  which  are 
long  with  respect  to  their  equivalent  flow  diameters  (fuselage  integrated  air  induction 
systems  typical  of  advanced  tactical  aircraft) .  For  relatively  long  air  induction  sys¬ 
tems,  it  becomes  necessary  to  divide  the  subsonic  duct  into  more  than  one  lumped 
volume.  Although  simulation  accuracy  increases  with  increasing  number  of  volume 
lumps  used,  practical  considerations  (computer  storage  and  run  time)  limit  the  num¬ 
ber  to  three  or  four.  All  the  models  reported  herein  divide  the  subsonic  duct  volume 
into  three  initially  equal  volumes . 


This  section  outlines  the  theory  and  equations  for  three  models,  all  based  on  the 
'lumped  parameter'  concept  for  the  subsonic  duct  downstream  of  the  NS.  Each  model 
treats  both  mixed  and  external  compression  air  induction  systems  at  supersonic 
speeds.  The  three  models  are: 

Model  1  —  Three  Lump  Standard  Model 

Model  2  —  Three  Lump  Buzz  Model 

Model  3  —  Three  Lump  Hammershock  Model 

The  2D  mixed  compression  (2DM)  air  induction  system  is  used  to  illustrate 
Model  1,  while  the  2D  external  compression  (2DE)  air  induction  system  is  used  to 
illustrate  Models  2  and  3 .  Since  only  one  air  induction  system  is  used  to  illustrate 
each  model  in  detail,  any  key  differences  in  the  basic  model  after  it  has  been  adapted 
to  the  other  air  induction  system  are  also  described.  Each  theoretical  model  i3 
divided  into  three  segments  —  Input  Data,  Initial  Conditions,  and  Dynamic  Calculations. 
The  Input  Data  segment  supplies  the  others  with  all  the  information  regarding  inlet 
steady  state  operating  conditions,  inlet  geometry  and  the  engine  inputs  to  the  air 
induction  system.  The  Initial  Conditions  segment  calculates  the  steady  state  air  induc¬ 
tion  system  performance  in  order  to  supply  the  Dynamic  segment  with  initial  values  of 
key  simulation  parameters  which  are  required  before  the  dynamic  calculations  can  be 
initiated.  The  complete  dynamic  problem,  including  the  response  of  the  relevant  con¬ 
trol  parameters  to  any  input  transient,  is  solved  in  the  Dynamic  segment.  Block 
diagrams  containing  the  structural  statements  required  for  each  model  are  included. 

A  detailed  discussion  of  the  lumped  parameter  concept  and  the  theoretical  approach 
used  in  formulating  each  model  follows . 

'LUMPED  PARAMETER'  CONCEPT 

The  application  of  this  concept  to  the  thermodynamic  system  can  be  explained 
with  the  help  of  an  analogous  electrical  network  with  lumped  parameters.  The  ther¬ 
modynamic  system  and  the  corresponding  electrical  analogy  are  shown  in  Figure  2. 
Inertia  (dynamic)  and  friction  (steady  state)  forces,  which  are  analogous  to  inductance 
and  resistance  in  the  electrical  system,  are  'lumped'  and  inserted  between  'lumped 
volumes'  with  stored  mass  and  energy,  the  latter  being  analogous  to  capacitance  in 
the  electrical  circuit. 
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Electrical  Analogy 


Input 

Voltage 


Figure  2.  Electrical  Analogy  for  ’Lumped  Parameter' 
Thermodynamic  System 
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The  thermodynamic  lumped  parameter  concept,  as  applied  to  supersonic  air 
induction  systems,  is  shown  in  Figure  3,  using  the  2 DM  inlet.  The  subsonic  duct  is 
first  divided  into  three  initially  equal  volumes,  Vol^.,  Vol^,  and  Vollg.  The  total 
pressures  and  temperatures  within  each  of  these  volumes,  during  the  simulation,  are 
calculated  using  the  equation  of  state  in  conjunction  with  the  unsteady  continuity  and 
energy  equations.  Detailed  derivation  erf  the  total  pressure  and  temperature  rate 
equations  is  presented  in  Appendix  I.  The  volume  erf  the  first  lump  (Volg)  changes 
with  the  movement  of  its  upstream  boundary  (represented  by  the  NS)  under  the  influ¬ 
ence  of  a  transient.  Boundaries  erf  the  remaining  two  volumes  (Vol.^  and  Vol^g)  are 
fixed.  The  mexlel  is  limited  to  cases  where  the  NS  does  not  cross  the  fixed  bound¬ 
aries,  since  this  situation  could  only  arise  when  the  bypass  control  malfunctions 
during  an  extremely  severe  transient,  a  rather  unlikely  occurrence  for  the  cases  of 
interest. 

To  calculate  rate  of  change  of  flow  rate  (Wg  and  W14)  at  each  of  the  two  fixed 
boundaries,  the  subsonic  duct  is  divided  into  two  halves,  each  half  containing  a  fixed 
boundary  (Stations  8  and  14) .  The  rate  of  change  of  flow  rates  at  these  boundaries  are 
then  determined  by  calculating  the  rate  of  change  of  momentum  of  each  volume  half 
due  to  the  instantaneous  net  imbalanced  force  acting  on  it  during  the  simulation.  For 
this  calculation  both  the  volumes  are  frozen  at  their  NS-at-throat  value  during  the 
transient.  The  single  fixed  boundary  for  this  calculation  is  located  at  the  volumetric 
center  of  Vol^  •  This  deriwtion  is  presented  in  Appendix  II. 

MODEL  1  -  THREE  LUMP  STANDARD  MODEL 

This  model,  as  applied  to  the  2DM  air  induction  system  (Figure  3),  is  limited  to 
simulating  its  response  in  the  started  mode.  Bypass  and  throat  area  controls  are 
provided  to  maintain  the  inlet  in  the  started  mode  and,  therefore,  for  transients 
severe  enough  to  cause  an  unscart  [excluding  buzz  (simulated  by  Model  2)  and  engine 
stall  (simulated  by  Model  3)],  the  control  system  response  may  have  to  be  modified 
to  avoid  the  unstart. 

This  model,  when  adapted  to  the  2DE  air  induction  system,  includes  an  upstream 
supersonic  model  for  simulating  NS  external  movement  on  the  ramp  upstream  of  the 
cowl.  A  description  of  this  upstream  supersonic  model  for  the  2DE  inlet  follows  the 


10 


detailed  discussion  of  the  theoretical  approach  in  formulating  the  2DM  air  induction 
system.  Model  1  is  capable  of  calculating  system  dynamic  response  to  both  external 
and  internal  transients.  However,  in  discussing  the  theoretical  approach,  an  internal 


transient  (rapid  changes  in  engine  flow  demand)  has  been  used  as  an  example. 

Dynamic  Calculations 

- r — — — - 

/ 

The  block  diagram  for  the  dynamic  calculations  is  presented  on  Sheets  1  to  4, 
Figure  4.  In  the  program,  the  NS  motion  (X)  during  a  transient  (external  or  internal), 
with  respect  to  the  duet,  is  calculated  using  the  equation 


*  (rt»  -»k) 


where 


c*  •  Jr  3 


=  the  Mach  number  of  the  flow,  with  respect  to  the  duct,  upstream  of  the 


and  Mx  =  the  Mach  number  of  the  flow,  with  respect  to  the  NS,  upstream  of  the  NS 


I  l±i  /£*y\ 

si  *  {rs* ) 


X  may  be  positive  or  negative,  depending  on  whether  the  NS  is  travelling  in  the  down¬ 
stream  or  upstream  direction,  respectively.  The  NS  is  assumed  to  be  without  inertia 
or  thickness.  Instantaneous  static  pressures  PgX  and  Pgy  across  the  NS  [Equa¬ 
tion  (2)1  are  calculated  independently  from  the  upstream  and  downstream  flow  proper¬ 
ties  as  follows: 

Properties  Upstream  of  the  NS.  The  instantaneous  value  of  Pgx  is  obtained  by 
first  calculating  the  flow  function  upstream  of  the  NS. 


A.  ,7tx 
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Figure  4  Continued:  Sheet  2  of  4 


Figure  4  Concluded:  Sheet  4  of  4 


where  is  the  local  duct  area  upstream  of  the  NS. 

is  then  obtained  using  this  value  of  flow  function  and  the  input  curve. 


Knowing  Tgv,  and  Pgx  can  then  be  calculated  as  follows. 


T*.  *  T„  /C  '  +  K  ) 


( r: ) 


Properties  Downstream  of  the  NS.  The  response  of  the  subsonic  duct  volume 
(and  hence,  the  instantaneous  value  of  PgY)  to  rapid  changes  in  flow  (Wg)  at  the  com¬ 
pressor  face  is  calculated  by  first  computing  the  rate  of  change  of  total  pressure  and 
temperature  within  each  of  the  three  volume  lumps,  Vol^g,  Vol12,  amI  Vo*g»  8tartta8 
from  the  compressor  face.  Thus,  within  Vollg,  the  total  pressure  and  temperature 
rate  equations  (see  Appendix  I  for  derivation)  can  be  written  as 
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Inflow  at  each  of  the  two  fixed  boundaries  is  computed  by  integrating  the  rate  of 
change  of  flow  (Wg  and  W^).  Thus,  at  the  fixed  boundary  station  14,  the  inflow  W^4 
into  Vol^  g  changes  at  the  rate 
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(9) 


■  rrc 


r  CpT,t  -  ^tz  -  °-saPu) 


where  I.,  is  the  dynamic  impedance  constant  (derived  in  Appendix  II) 

««*  Af>_  =  k„  <.PTX  -  P**) 

=  o  oz6(,Prx  -  P<»x) 


Kj^  is  an  empirical  subsonic  diffuser  pressure  loss  factor  which  is  a  function  of  M^. 
The  value  of  .026  selected  for  this  factor  by  Reference  2  in  their  simulation  models  is 
also  used  here. 

Similarly,  within  Vol19,  ^  TTl2  are  °btalnec*  from  equations 


rTI2-  = 


12.  K 

-  ( 
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"pTii.  *^<<4  ) 
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where,  at  the  fixed  boundary  station  8,  the  inflow  Wg  into  Vol12  changes  at  the  rate 


L  (?TC  - 


T\Z  - 


The  total  pressure  and  temperature  rate  equations  for  Volg  include  the  effect  of 
the  moving  upstream  boundary,  represented  by  the  NS,  so  that 

Ko  =  JLf: _ [ttxO'x  -  wT6y  +  wxSK)_Tr&w(.l  *  il  <“>) 

l  J 


Tn  pi  fepi.  («&  *  wT6y  -  <15 

Ptcp  L  VL  vu  j 

Inflow  into  Volg,  where  the  upstream  boundary  rr^resented  by  the  NS  is  moving,  can 
be  visualized  to  consist  of  two  parts: 

(a)  Inflow  with  respect  to  the  duct 

(b)  Inflow  swept  by  the  upstream  boundary  area  A^,  moving  at  velocity  X  rela- 
tive  to  the  duct,  through  air  at  density  p^,  which  is  equal  to  Pg^  A^.  X. 
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Thus,  inflow  with  respect  to  the  normal  shock 

-  A*  * 

=■  4*  l^XSK 

where,  by  definition, 


Ul 


XSK 


Ax  X 
*  Tsx 


and  is  negative  for  positive  X  (downstream  NS  movement)  and  vice-versa. 


(16) 


The  total -to-statlc  pressure  and  temperature  ratio  terms  C'  and  Z  (defined  in 
Appendix  I)  are  obtained  by  entering  the  input  data  tables  C:g  and  Zg,  repeatedly  with 
the  following  flow  functions  based  on  the  instantaneous  mean  flow  rates  within  each 
volume  lump. 

(17) 

(18) 

(19) 

Using  PT6  (obtained  by  integrating  Equation  14),  Pgy  is  next  calculated  as  follows. 

The  flow  function  Fy  may  be  written  as 


C-  ~ 
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From  this  ratio,  Pgy  is  given  by 


The  instantaneous  values  of  PgX  and  Pg^  (obtained  from  Equations  6  and  22,  respec¬ 
tively)  are  then  used  in  Equation  (2)  to  obtain  for  calculating  the  normal  shock  (NS) 
velocity  X  (Equation  (1)]. 

Air  Induction  System  Controls:  Control  logic  is  provided  in  this  model  to  pre¬ 
vent  two  types  of  unstarts  associated  with  mixed  compression  inlets:  Choked  throat 
unstart  and  the  unstart  caused  by  disgorgement  of  the  NS  resulting  from  reduction  in 
engine  demand.  Choked  throat  unstart  occurs  when  the  throat  cannot  pass  the  flow 
being  supplied  by  the  inlet,  resulting  in  an  unstable  normal  shock  appearing  just 
upstream  of  the  throat  (in  addition  to  the  NS  located  downstream)  which  is  immediately 
disgorged  to  allow  subsonic  spillage  of  the  excess  air.  Choked  throat  unstart  can  be 
caused  by  (a)  increase  in  the  throat  flow  function  towards  the  choked  value  due  to 
reduction  in  MQ,  (b)  reduction  in  throat  area,  and  (c)  increase  in  the  flow  being  sup¬ 
plied  to  the  inlet  due  to  increase  in  angle  of  attack.  Thus,  during  this  type  of  an 
unstart,  two  normal  shocks  can  instantaneously  exist  in  the  induction  system,  although 
the  choked  throat  unstart  will  probably  also  cause  the  supercritical  NS  (i.e. ,  the  ter¬ 
minal  shock)  to  move  rapidly  upstream,  overtake,  and  ultimately  coalesce  with  the 
unstable  normal  shock  originating  at  the  throat  before  the  latter  reaches  the  cowl  lip. 
In  the  program,  a  variable  area  bypass  is  used  to  control  the  terminal  shock  motion 
and  the  inlet  throat  area  is  programmed  to  open  if  the  throat  tends  to  choke.  The  con¬ 
trol  logic  flow  is  presented  in  Figure  4,  Sheet  4. 

Bypass  Area  Control:  The  basic  bypass  area  control  error  function  con¬ 
sists  of  the  following  four  error  signals: 

^Dx  •+■  ^Ae*  Ae>y  (23) 

where,  C^x  AX  is  the  proportional  signal  on  shock  position  deviation 

"  <*REF  -  *> 

GDx  “X  is  the  derivative  signal  on  shock  position  deviation,  AX 
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nd  ANp  is  the  derivative  signal  on  engine  fan  speed  deviation 

ANp  =  (NpREF  -  Np).  Since  the  engine  fan  speed 
directly  effects  the  engine  flow  demand, this  is  essentially 
an  engine  flow  (W2E)  derivative  signal. 

G^b  ACy  is  the  derivative  signal  on  bypass  area,  Agy,  and 

GAX*  GDX*  GDN  an<1  GAB  are  the  gains  ussociated  with  the  signals  defined 
above. 


Also, 


is  the  derivative  lag  on  shock  position  deviation,  AX, 
is  the  derivative  lag  on  engine  fan  speed  deviation,  ANp.and 

is  the  derivative  lag  on  bypass  area,  Agy 


Each  of  the  above  four  error  signals  provide  an  overall  bypass  area  error  signal 
which  is  a  function  not  only  of  the  position  error  (XREp  -  X),  but  also  of  the  rate 
at  which  the  deviation  of  X  and  NR  occurs  from  their  respective  reference  val¬ 
ues,  and  the  rate  at  which  Agy  opens  or  closes. 

A  derivative  lag  is  a  derivative  of  the  lagged  input  and  is  obtained  by  subtracting 
the  lagged  input  from  the  input  and  dividing  the  result  by  the  time  constant  r. 
Thus,  using  Laplace  transforms 


2 

i 
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£  S 
1  +  fS 


(24) 


The  expression  on  the  RHS  in  the  above  equation  is  the  derivative  of  the  lagged 
2 

input  ~+Tg  as  defined  by  the  theorem  on  Laplace  transforms  of  derivatives  (Ref¬ 
erence  7). 


Each  of  the  derivative  signals  in  the  error  function  (23)  is  calculated  in  the 
model  as  follows,  assuming  a  practical  value,  in  terms  of  controls  hardware 
response,  of .  01  for  the  time  constant. 

The  derivative  signal  on  shock  position  deviation  (AX)  is 


Gox  Ax 


(AX 


A  x 

l  -+•  oi  s 


(25) 
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Note  that  AX  is  negative  for  downstream  movement  of  the  normal  shock.  This 
essentially  means  that  the  signal  will  tend  to  increase  the  negative  portion  of  the 
total  error  in  the  bypass  logic  if  the  rate  of  deviation  increases  as  the  NS  moves 
downstream,  thus  tending  to  close  the  bypass  at  a  faster  rate  to  reduce  the 
deivation. 

Similarly,  the  derivative  signal  on  engine  fan  speed  deviation  (ANp)  is 


=  G 


•Dn J 


fO  p  \  | 


(26) 


If  the  fan  speed  increases  beyond  Np^pp,  the  engine  instantaneously  demands 
more  flow  and  the  NS  again  tends  to  move  downstream.  With  ANp  negative,  an 
error  signal  similar  to  AX  is  input  to  drive  the  NS  back  towards  the  throat. 


The  AgY  derivative  signal  is 
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(27) 


This  is  a  damping  input  on  the  closing/opening  rate  of  AgY  itself  to  avoid  large 
overshoots.  The  derivative  input  (i.e. ,  the  expression  in  brackets  in  equation 
(27)]  is  deliberately  reversed  to  allow  for  a  positive  error  when  ABY  closes, 
(thus  damping  the  overall  negative  error)  and, correspondingly,  a  negative  error 
when  it  opens. 


Throat  Area  Control:  The  throat  area  control  logic  sets  up  an  error  func¬ 
tion  (A^g^^  -  A^,)  for  any  deviation  of  the  actual  throat  area  (A^,)  from  the 
scheduled  value  (A^g^)  and  tries  to  hold  the  scheduled  value  by  reducing  the 
error.  Thus,  in  trying  to  hold  the  scheduled  value,  the  control  not  only  accom¬ 
plishes  the  task  of  avoiding  a  choked  throat  unstart  but  also  tries  to  maintain 
inlet  performance  as  close  to  optimum  as  possible  during  a  transient. 


Initial  Conditions 


This  segment  supplies  the  dynamic  segment  with  the  initial  values  of  the  follow¬ 
ing  parameters,  calculated  from  the  input  data,  which  are  required  before  the  dynamic 


calculations  can  be  initiated:  (X)IC,  (PTg)lC,  (PTl2>JC.  (PT18>lC,  (T' 


T6^1C’ 


Vl2*lC’  ^TT18^IC’^Wq^ 


8'IC’  ^14*10’  ^ 


by>ic»  ('Vic* 
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This  is  achieved  by  calculating  the  air  induction  system  steady  state  performance  at 
the  given  operating  condition  as  shown  in  the  block  diagram  (Figure  5).  The  NS  initial 
condition  (X)j^  is  obtained  from  the  free  stream  Mach  number  using  the  input  data 
table  =  f(MQ).  The  initial  values  of  various  flow  rates  !(Wq)j£.»  (W^)^., 

("tBY^C1,  areas  HAxhc’  ^T^IC*  311(1  the  NS  uPstream  total  temperature  fI’T^)IC 
are  calculated  next, using  the  input  data  tables  as  shown  in  Sheet  1.  These  are  then 
used  as  inputs  in  Sheet  2  to  calculate  properties  behind  the  NS.  (P-pg)^  is  obtained 
from  the  total  pressure  drop  corresponding  to  the  NS  upstream  Mach  number, 

(Mx)jc*  The  initial  total  pressures  (P-p^IC  anc*  ^PTi8^IC  are  tlien  ol)tained  as 
follows. 


(*T  12  )ic  c  )  ic  "  c>  •  S'  (  A  Pu  )  Jc 

T»<i)tc  =.  (PtOic  =  CPT^)ic  "  £  C 6  )  I  c 
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where 
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(APLVc  «  C'OL  **  [  c  Prx  )lc  " 

The  total  temperatures  throughout  the  air  induction  system  are  initially  constant  and 
equal  to  the  initial  freestream  total  temperature  (TTQ  )IC  since  the  flow  process  is 
assumed  to  be  adiabatic. 

At  the  initial  condition,  the  bypass  area  is  sized  such  that  the  NS  is  located  at 
the  desired  supercritical  location  as  follows-. 
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o.  s  <h  (  P n  )  ic 
The  bypass  area  is  assumed  to  be  choked  during  the  entire  simulation  run. 

Input  Data 

Input  data  required  for  the  2DM  air  induction  system  dynamic  s  mulation  is  pre¬ 
sented  in  Table  I.  The  functional  relations  are  shown  for  each  input  variable  along 
with  the  figure  in  which  it  is  presented.  The  overall  data  is  input  into  the  program 
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igure  5.  Model  1  -  2DM  Three  Lump  Standard  Model  Block  Diagram 
Initial  Conditions:  Sheet  1  of  2 


ShetZl 


TABLE  I.  MODEL  1  -  2DM  THREE  LUMP  STANDARD  MODEL  INPUT  DATA 


INPUT  PARAMETER 

FUNCTIONAL  RELATION 

FIGURE  NO. 

Inlet  Key  Operating  and 
Geometric  Parameters 

* 

0 

i'O  or  constant 

P 

constant 

- 

* 

H 

f(t)  or  constant 

- 

★ 

uo 

f(t)  or  constant 

- 

acap 

constant 

- 

P0STD 

f(H) 

- 

T0STD 

f(H) 

- 

ml 

m0>  «> 

13(a) 

mr 

f<ML*  aR2) 

13(b) 

PTL/PT0 

W0>  «) 

13(c) 

PTR/PTL 

f^L’  °W 

13(d) 

ao/acap 

f(MQ,  a) 

13(e) 

°R2 

£(M0) 

13(f) 

Ax 

13  (k) 

atsch 

f<QR2> 

13(g) 

wbr/wo 

f(ML) 

130 

WHX/W0 

mL>  x) 

130 

WBi-/W0 

mv  x) 

130 

XSU 

f<QR2) 

130 

XSL 

^°R2) 

130 

XREF 

^RZ* 

13(d) 

psy/Pty 

f(Fy) 

13  6>) 

MXI 

fOPja'  yfP-t 

13  (n) 

Scale 

Constant 

- 

i 


f 


TABLE  I.  MODEL  1  -  2DM  THREE  LUMP  STANDARD  MODEL  INPUT  DATA 

(Continued) 


INPUT  PARAMETER 

FUNCTIONAL  RELATION 

FIGURE  NO. 

Lumped  Volume  Parameters 

(PT6/PTX)IC 

mx) 

13(o) 

A12 

constant 

- 

A18 

constant 

- 

X8 

constant 

- 

v°i6 

i(oR2,  X) 

130) 

v°i12 

constant 

- 

oo 

r-4 

o 

> 

constant 

- 

!1 

constant 

- 

*2 

constant 

- 

V  Z12’  Z18 

f(Fg)»  f<F12),  f(P18) 

13  (r) 

CV  C*12*  C *  18 

f(F6),  f(F12),  f(Flg) 

13  (q) 

. 

Engine  Parameters 

* 

nf 

f  (t)  or  constant 

- 

* 

W 

W2EC 

f^p/'/fig)  or  constant 

13  (m) 

*These  parameters  define  specific  input  transients.  One  or  more  of  these  is  input 
as  a  function  of  time  for  a  particular  simulation  run,  the  others  being  constant. 


! 
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either  as  a  constant  or  as  one  or  two  dimensional  tabulated  data.  Parameters  marked 
with  an  asterisk  define  specific  input  transients.  Any  one  or  more  of  these  may  be 
input  as  a  function  of  time  for  a  particular  simulation  run,  the  others  remaining  con* 
stant.  This  concludes  the  discussion  on  Model  1  as  applied  to  the  2DM  air  induction 
sysiem.  The  key  difference  in  this  model,  when  it  is  used  for  the  2DE  air  induction 
system,  is  the  model  for  simulating  NS  travel  on  the  inlet  ramp.  This  is  discussed 
below. 


Standard  Supersonic/ Subsonic  Upstream  Model 

For  the  2DE  air  induction  system,  a  model  is  required  to  simulate  NS  travel  on 
the  ramp  when  the  inlet  is  operating  subcritically.  The  approach  is  shown  schemat¬ 
ically  in  Figure  6.  Figure  6(a)  shows  the  NS  within  the  cowl  when  the  inlet  operates 
supercritically  at  an  off-design  Mach  number  and,  hence,  with  only  supersonic 
spillage  (W^p  -  Wq).  The  inflows  and  outflows  from  the  first  lumoed  control  vol¬ 
ume,  Volg,  when  the'  inlet  operates  subcritically,  are  shown  in  Figure  6(b).  The  sub¬ 
sonic  spillage  Wgp  required  to  calculate  the  flow  entering  the  inlet  (WQ  -  Wgp)  i8 
obtained  from  Moeckel's  continuity  model.  This  model  defines  the  relationship 
between  shock  stand-off  distance  (L)  and  the  subsonic  spillage  and  is  explained  in 
detail  in  Appendix  III. 


MODEL  2  -  THREE  LUMP  BUZZ  MODEL 

A  number  of  theories  are  available  in  literature  which  attempt  to  explain  the 
rather  complex  unsteady  phenomena  called  inlet  "buzz"  (References  8  and  9).  Inlet 
buzz  generally  occurs  at  subcritical  operation  when  the  inlet  mass  flow  ratio  falls 
below  a  certain  level  causing  large  excursions  of  the  normal  shock  (NS)  fore  and  aft 
along  the  inlet  compression  surfaces  at  a  high  frequency.  The  resulting  large  periodic 
total  temperature /pressure  variations,  and  the  associated  dynamic  distortion  at  the 
compressor  face,  can  lead  to  engine  stall. 

The  analytical  model  described  herein  does  not  attempt  to  predict  conditions  at 
which  buzz  occurs.  Rather,  its  function  is  to  indicate  the  buzz  intensity  (i.  e.,  frequency, 
NS  excursion  distance,  PT2  and  TT2  amplitudes,  etc.),  once  it  has  been  triggered  at 
a  particular  inlet  operating  condition.  The  buzz  concept  adopted  for  this  model,  and 
the  theoretical  approach  used  in  simulating  it,  is  discussed  below  using  the  2DE  air 
induction  system.  The  same  model  applies  to  the  2DM  air  induction  system  except 
for  the  added  internal  NS  travel  from  the  throat  to  the  cowl  defined  as  the  'unstart ' 
portion  of  the  transient. 
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Capture  Flow  Area 
Critical  Flow  Area 


Figure  6.  Model  1  -  Standard  Supersonic/ Subsonic  Upstream  Model 
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The  Buzz  Concept 


The  theory  for  burr  adopted  here  focuses  on  the  massive  boundary  layer  separa¬ 
tion  due  to  the  normal  shock-boundary  layer  interaction  on  the  inlet  compression 
surfaces  as  the  source  of  instability  (Reference  2).  The  criteria  for  normal  shock 
induced  boundary  layer  separation  is  presented  in  Reference  10.  It  states  that  separa¬ 
tion  occurs  when  the  static  pressure  across  the  NS  rises  to  approximately  1.89.  This 
static  pressure  rise  requires  that  the  Mach  number  entering  the  shock  be  1.33  or  above. 

The  buzz  concept  is  shown  schematically  in  Figure  7.  The  turbulent  separated 
boundary  layer  region  contains  relatively  low  energy  air  compared  to  that  outside  the 
boundary  layer.  This  low  energy  turbulent  separated  boundary  layer  zone  tends  to  act 
like  a  wedge  superimposed  on  the  inlet  compression  surface,  thus  throttling  the  flow 
entering  the  inlet  and  causing  significant  quantities  of  the  high  energy  air  to  spill  sub- 
sonically  outside  the  cowl. 

Figure  7(a)  schematically  shows  first  lumped  control  volume  (Volg),  and  the 

inflow  and  outflow  through  the  boundaries.  The  upstream  boundary  is  the  NS  which 

extends  up  to  the  critical  flow  streamline.  The  inflow  (W^)  is  compressed  somewhat 

by  the  separated  boundary  layer  as  it  enters  the  NS  through  the  effective  area  A£X, 

which  is  the  geometric  area  (A^)  less  the  separated  boundary  layer  area  (AA£X). 

This  inflow  includes  the  subsonic  spillage  flow  (Wgp).  The  control  volume  includes 

both  the  high  and  the  low  energy  streams,  with  the  latter  merely  acting  as  a  throttling 

agent,  allowing  only  the  throttled  throat  flow  W  to  exit  (or  enter,  if  reverse  flow 

ET 

conditions  exist)  through  the  fixed  boundary  at  Station  8.  The  remaining  outflows  from 
the  volume  are  the  ramp  and  throat  bleeds  and  the  subsonic  spillage  between  the  NS 
and  the  cowl. 

The  key  parameters  selected  for  simulating  inlet  buzz  are  the  effective  throat 
area  ratio  A£T/A^,  and  the  total  pressure  ratio  (Figure  7(b)  and  (c)).iAs  the  NS 

travels  upstream,  separation  increases  significantly,  causing  massive  throttling  when 
the  NS  reaches  the  first  ramp.  Figure  7(b)  shows  that  the  effective  throat  area  ratio, 
AET/Ap,  decreases  to  0.3  for  the  2DE  inlet  operating  at  Mach  2.5,  due  to  throat 
throttling  as  the  NS  travels  upstream  from  the  cowl.  The  decrease  in  Aet/At  is 
based  on  throat  blockage  levels  used  in  Reference  11.  which  comoares  XB-70 
flight  test  data  with  simulation  predictions  of  inlet  unstart  and  buzz. 

The  throttling  of  flow  due  to  massive  boundary  layer  separation  during  the  buzz 
cycle  causes  significant  instantaneous  flow  imbalance  within  the  subsonic  duct.  That 
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is,  during  the  separated  flow  portion  of  the  buzz  cycle,  inflow  to  the  subsonic  duct  is 
drastically  reduced  while  the  engine  demand  remains  unchanged,  resulting  in  a 
sharp  reduction  in  pressure  throughout  the  duct.  As  the  pressure  drops  below  a 
critical  value  (specifically  when  n/P~,„  -  .315),  reattachment  of  the  flow  occurs, 

x  iZ  XU 

thus  increasing  the  throat  flow  area  and  the  subsonic  duct  pressure  level.  As  the 
pressure  rises  above  a  critical  level  (specifically,  when  P>pi2^PT0  sePara" 

tion  of  the  boundary  layer  is  reinitiated.  The  reattachment/separation  critical  values 
of  .315  and  .519,  respectively,  were  established  in  Reference  2. 


Dynamic  Calculations 

The  overall  calculation  procedure  for  the  buzz  condition  is  the  same  as  for  Model 
1.  That  is,  the  NS  velocity  X  is  obtained  from  the  instantaneous  static  pressures  PgX 
and  PgY  calculated  independently  from  the  flow  properties  upstream  and  downstream  of 
the  NS.  However,  this  model  differs  in  the  method  used  to  calculate 

(a)  properties  upstream  of  the  NS,  and 

(b)  properties  within  the  first  lumped  volume  (Volg). 

Also,  additional  logic  has  been  added  for 

(a)  simulating  buzz  using  a  boundary  layer  separation/reattachment  model, 
with  on/off  and  hysteresis  characteristics,  that  is  controlled  by  key  buzz 
parameters,  and 

(b)  calculating  dynamic  conditions  downstream  of  the  NS  for  buzz  severe  enough 
to  cause  flow  reversal  within  the  subsonic  diffuser. 

A  block  diagram  of  the  equations  used  in  formulating  the  buzz  model  is  presented 
in  Figure  8,  Sheets  1  to  6.  The  key  differences  and  additions  mentioned  above  are 
next  discussed  with  the  help  of  this  block  diagram. 


Properties  Upstream  of  the  NS.  The  instantaneous  static  pressure  P  is 
’  „  sx 

obtained  by  first  calculating  the  appropriate  NS  upstream  Mach  number  (Mx)  by 

entering  the  input  data  curve 


M*  ;  f(F„) 

with  the  corresponding  value  of  the  flow  function  (Fv)  defined  as 

X 


‘'X 


U)  jc  4  X. 


(32) 


(33) 
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2DE  Three  Lump  Buzz  Model  Block  Diagram 
Dynamic  Section:  Sheet  1  of  6 
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Figure  8  Continued:  Sheet  2  of  6 
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Continued 
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Figure  8  Continued:  Sheet  4  of  6 
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Figure  8  Continued:  Sheet  5  of  6 


Figure  8  Concluded:  Sheet  6  of  6 


The  effective  flow  area  A„„,  upstream  of  the  NS,  is  obtained  by  entering  the  input 
data  curve 

Ae*  =  f  C*)  (34) 

with  the  corresponding  instantaneous  value  of  X.  This  curve  defines  the  effective 
flow  area  upstream  of  the  NS  for  subcritical  inlet  operation  (NS  outside  the  cowl, 

X  <  Xc)  and  the  duct  cross-sectional  area  upstream  of  the  NS  for  critical/super¬ 
critical  inlet  operation  when  the  NS  is  inside  the  cowl  (X  >  Xc). 

For  adiabatic  flow  the  total  temperature  upstream  of  the  NS  (T^,^)  is  equal  to 
the  freestream  total  temperature  T^^.  This  temperature  and  are  used  to  calculate 
the  corresponding  static  temperature 


Tsx  =  TYx/(l  + 

PgX  is  then  calculated  from  the  relationship 


*  *Tc 


where  the  total  pressure  ratio  p-j>X^T0  is  ky  entering  the  input  data  curve 

PTX/PTo  =  i(x)<  with  the  corresponding  instantaneous  value  of  X. 

Properties  Downstream  of  the  NS.  Before  discussing  the  response  of  the  sub¬ 
sonic  duct  volume  to  the  buzz  transient,  it  is  appropriate  to  examine  the  throat  throttl¬ 
ing  logic  that  simulates  the  buzz  transient  (Figure  8,  Sheet  3).  In  calculating  the 
throttled  throat  flow,  it  is  assumed  that  the  instantaneous  flow  (Wg)  exiting  (for  forward 
flow)  or  entering  (for  reverse  flow)  Volg  at  Station  8  is  equal  to  the  instantaneous 
throttled  throat  flow,  W£T.  Secondly,  the  lumped  properties  at  Station  8  are  used  in 
calculating  Wg.  This  is  permissible  since  the  'lumped  parameter'  concept  assumes  that 
properties  within  the  lumped  volume  Volg)  and  hence,  the  throat  properties,  are 
instantaneously  equal  co  the  property  values  at  the  exit  Station  8.  The  rate  of  change 
of  throttled  throat  flow  (Wg)  is  calculated  next  as  follows. 


Calculation  of  W^.  The  throat  blockage,  (AA^), is  first  calculated  by  enter¬ 
ing  the  input  data  table 

A  =  -f  (  x )  (37) 
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with  the  corresponding  instantaneous  value  of  X.  Then, 


Now,  referring  to  Figure  7(c),  the  flow  once  separated  remains  separated  as 
long  as  P-p^/^TO  *s  a'aove  and,  therefore,  for  separated  flow  the  full  value 
of  the  blockage  (AA£T)  is  used 


A  Aer 


I .  C  A  A  ET 


When  P'pi2^T0  ^roPs  below  .315,  attachment  occurs  and  flow  remains  attached 
until  the  pressure  ratio  rises  above  .510.  During  this  period  when  flow  is 
attached,  blockage  is  assumed  to  be  two  percent  of  its  separated  flow  value,  so 
that  for  attached  flow 

A  l\Er'  --  c.oi.AAET  (40) 

Hence,  the  effective  throat  area,  at  any  instant,  can  be  calculated  by  using  the 
appropriate  AA^  for  either  separated  or  attached  flow  in  the  equation 

Ae_  .  AT  -  AAe't  (41) 


The  throat  throttling  flow  functions  for  both  forward  flow  (FTFWD)  and 
reverse  flow  (FTREy)  ar8  next  ca^cu^ed  by  entering  the  input  data  curve 
[Ft  =  f(PT/Pg)  ]  with  the  corresponding  value  of  the  appropriate  total-to-static 
pressure  ratio  at  Station  8.  The  procedure  is  discussed  below. 

For  forward  flow,  the  total-to-static  pressure  ratio  in  Volg  (PTg/pgg)  is  used  to 
calculate  as  follows. 

The  input  data  curve 

=  -f  (F12.)  (42) 

Pt.r- 

includes  the  front  half  of  the  duct  steady  state  friction  losses  within  the  steady 
state  pressure  recovery  term,  DurinS  simulation,  the  above  input 

curve  is  entered  with  the  corresponding  mean  flow  function  for  Vol^, 

fa  -  o  s  C  JTrt2-  (43) 

A  u-  Frri- 
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to  obtain  the  instantaneous  P^/P^,^  value.  This  ratio  and  the  instantaneous 
total  pressures  PTg  and  are  use<*  the  instantaneous  total-to- 

static  pressure  ratio 


•  (&Xfe) 


Then  the  input  data  curve  [FT  -  f(PT/Pg  ]  is  entered  with  the  pressure  ratio 


PT6/ps6  t0  obtain 


FrpuD  -  < 

Similarly,  for  reverse  flow,  the  same  curve  [Ft  =  f(P,p/Pg  ]  is  entered 
again  with  the  pressure  ratio  p-ri2^S12  to  °kta*n 

^rR£v  *  ^  (  pp—V  * 

\  Pit u/ 


Then,  for  forward  flow,  the  throttled  throat  flow  is 

U5  =  (^TFUiO  Ae.T  )/J  T-rfc  if  (47) 

and  for  reverse  flow,  it  is 

U»  -  (Fraev  ^  **r)/J*n  “  (48) 

Finally,  Wg  is  obtained  by  imposing  a  lag  on  the  appropriate  Wg,  using  a  typical 
value  of  .  005  for  the  time  constant  r^g. 

U«  -  (We  -  W*)/*w»  <«> 

.  K  -  *0/  •CoS"  (50) 

Calculation  of  the  Instantaneous  P^,.  The  response  of  the  subsonic  lumped 
volume  Volg  (and  hence,  the  instantaneous  value  to  PgY)  to  the  throttling  of  the 
flow  at  the  fixed  boundary  exit  Station  8  is  calculated  using  the  PTg  and  TTg 
rate  equations. 

For  forward  flow,  Wg  is  greater  than  zero  and  exits  Volg  at  the  instantaneous 
total  temperature  TTg,  so  that 

|JU>  *  [Trx  (w *  +  <)  -  rTfc  (u;t e>y  **•  usp  +  we)J  +  P— ^  -  (51) 


p  _  12.  R 

i  rfc>  _ 
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For  reverse  flow,  Wc  is  less  than  zero  and  enters  Vol„  at  the  instantaneous 
8  6 

total  temperature  TT12,  so  that 

The  appropriate  value  of  PTg  corresponding  to  forward  or  reverse  throttled  flow 
direction  (Wg  Is  positive  or  negative,  respectively)  Is  selected  through  program 
logic  (see  Figure  8,  Sheet  2). 

For  both  forward  and  reverse  flow, 

+  WTay  -  *****  1(58) 

PTtc.  L  Vk,  vg*  J 

Note  that  the  inflow  (W^)  through  the  NS  (Volg  upstream  boundary)  includes  the 
subsonic  flow  spillage  (Wgp)  when  the  NS  is  upstream  of  the  cowl,  Moeckel's 
model  (Appendix  III)  has  been  incorporated  into  the  buzz  model  to  calculate  this 
spillage  flow  during  the  simulation.  The  NS  standoff  distance  L  is  known,  at  any 
instant,  during  the  dynamic  calculations 


a  Xc  -  X 

foe 

X 

< 

(54) 

*  o 

for 

X 

* 

The  spillage  coefficient,  Kgp,  corresponding  to  the  NS  upstream  Mach  number 
M^,  is  obtained  from  Moeckel's  curve 

^  -  f  (M*).  (55) 

From  this,  the  corresponding  spillage  is 


U0  l 

*<sp 


(56) 


where  WQ  is  the  critical  inlet  flow. 

Using  PT6  (obtained  from  Equation  51  or  52)  PgY,  and  finally  X  are  calcu¬ 
lated  as  already  explained  for  Model  1.  The  detailed  procedure  is  shown  in 
Figure  8,  Sheet  2. 


Voijn  and  Vol^  Instantaneous  Total  Conditions.  The  instantaneous  total 
pressures  and  temperatures  within  the  volumes  Vol10  and  Vol.Q  are  obtained  by 

L&  lo 

numerically  integrating  the  corresponding  rate  equations.  In  Vol10,  since  for- 
ward  or  reverse  flows  can  occur  at  both  the  fixed  boundary  Stations  8  and  14, 
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any  one  of  four  combinations  of  forward/reverse  flows  can  exist  at  any  instant. 
Matching  of  these  flows  with  the  appropriate  instantaneous  total  temperatures, 
at  each  fixed  boundary,  leads  to  the  following  four  total  pressure  rate  equations: 


If  Wg  and  W]4  >  0, 


12.  * 


C  '  ?•  V 
iu  vu. 

If  Wg  >0  and  WJ4<  0 


tt<c  ~  w  it* 


] 


(57) 
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12. 
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(58) 


(59) 


(60) 


The  appropriate  value  of  Pmo  corresponding  to  the  instantaneous  WQ  and  W, 
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flow  directions  is  selected  through  program  logic  (see  Figure  8,  Sheet  4). 


For  both  forward  and  reverse  flows, 

T, 


T 


Tit 


_  lTl2. 

Ptiv 


,  -  **)  +  pT,v.l 


(61) 


The  calculation  procedure  for  obtaining  the  total  pressure  and  temperature 
rates  within  Vol-g,  for  forward  and  reverse  flows  at  the  two  boundary  stations 
14  and  2,  is  similar  to  that  discussed  above  for  Voi12  and  is  shown  in  Figure  8, 
Sheet  5.  An  option  is  provided  in  program  logic  for  calculating  FT18  with  reverse 
flow  at  the  compressor  face  (W2  <  0;  i.  e.,  during  severe  engine  stall).  An  appro¬ 
priate  matching  value  of  engine  stall  total  temperature  (TTES)  is  required  as  in¬ 
put  for  simulating  such  a  transient. 


Calculation  of  W^.  The  forward  (W^p^)  and  reverse  (W14REy)  rates 
of  change  of  flow  rates  at  the  second  fixed  boundary  (Station  14)  are  calculated 
by  applying  the  momentum  equation  to  the  half  duct  volume  containing  the  second 
fixed  boundary. 
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The  two  equations  are 

W«FW>  =  Y  [>W  -  ~  <="S  (~;  )  PT'»-]  (62) 

W,U«V  -  ~~  [?TU-  -  Pt‘6  *■  C '*  (-^)  fTls  J  (63) 

where  APL/PT  is  the  subsonic  diffuser  total  pressure  loss  due  to  friction,  and 
is  assumed  to  be  constant. 


Note  that  for  reverse  flow,  the  frictional  force  acts  in  the  upstream  direction. 
Hence,  the  change  in  sign  for  the  third  expression  in  brackets  in  Equation  (63). 


Thus, 


u 

=  WI4FU£) 

W,4  > 

*  u  i4iiev 

■V 

W,4  <  C 

(64) 


Initial  Conditions 

The  initial  values  of  thq  following  parameters  are  supplied  by  this  segment  to 
the  dynamic  segment  before  dynamic  calculations  are  initiated.  (X)^,  (P^)^, 

^T^IC’  <TTehc  (TT12)IC’  *TT18*IC*  ^IC’  ^14*10  * 

The  calculations  are  shown  in  Figure  9,  and  are  self-explanatory.  With  no  by¬ 
pass,  all  the  excess  flow  is  spilled  subsonically  at  the  cowl  lip  with  the  NS  located  at 
the  proper  initial  stand-off  distance  from  the  cowl  to  allow  for  the  required  spillage. 

Input  Data 

The  key  input  parameters  required  for  the  2DE  air  induction  system  buzz  simu¬ 
lation  are  presented  in  Table  II.  This  model  does  not  require  any  input  transients. 
Buzz  is  triggered  through  internal  program  logic  as  explained  in  the  dynamic  segment. 


MODEL  3  -  THREE  LUMP  HAMMERSHOCK  MODEL 

Rapid  cutoff  of  engine  demand,  usually  due  to  compressor  stall,  causes  an 
extreme  instantaneous  flow  imbalance  within  the  diffuser.  With  the  diffuser  outflow 
drastically  reduced  and  the  inflow  unaffected,  the  instantaneous  excess  of  inflow  over 
outflow  causes  a  sharp  increase  in  pressure  near  the  compressor  face.  With  the  dif¬ 
fuser  exit  blocked,  this  local  pressure  wave  increases  in  strength  and  propagates 
upstream  as  a  shock  front  or  hammershock  (HS).  The  pressure  buildup  within  the  duct 
is  relieved  through  subsonic  spillage  when  the  hammershock  exits  from  the  duct  at  the 


TABLE  H.  MODEL  2  -  2DE  THREE  LUMP  BUZZ  MODEL  INPUT  DATA 


INPUT  PARAMETER 

FUNCTIONAL  RELATION 

FIGURE  NO. 

Inlet  Key  Operating  and 
Geometric  Parameters 

a 

constant 

constant 

- 

H 

constant 

- 

M0 

constant 

- 

acap 

constant 

- 

ao/acap 

constant 

- 

°R2 

constant 

- 

VW2 

constant 

- 

APL/p 

constant 

- 

(Msp)ic 

constant 

- 

P0STD 

f(H) 

- 

T0STD 

f(H> 

- 

PTX/PT0 

f(X) 

22(a) 

^etMt 

f(X) 

22(d) 

Ax 

f(X) 

22(c) 

aex 

m 

22(c) 

KSP 

{(Mxe) 

16(e) 

Mxi 

f(Fx*  /R) 

13(n) 

PSY/PTY 

f(Fy) 

13<p) 

Scale 

constant 

46 


TABLE  II.  MODEL  2  -  2DE  THREE  LUMP  BUZZ  MODEL  INPUT  DATA 

(Continued) 


INPUT  PARAMETER 

FUNCTIONAL  RELATION 

—  . -  1  " 

FIGURE  NO. 

Lumped  Volume  Parameters 

(PT6/pTx>ic 

f(Mx) 

13(o) 

A12 

constant 

- 

A18 

constant 

- 

X8 

constant 

- 

V0l6 

f(X) 

22(b) 

Voll2 

constant 

- 

V0l18 

constant 

- 

1 2 

constant 

- 

PS6//PT12 

^12* 

16(f) 

PT 

f(PT6/PS6) 

16(g) 

Z6’  Z12’  Z18 

f(F6),  f(F12),  f(F18) 

13  (r) 

C6'  C'!2'  C’l8 

f(F6),  f(F12),  f(F18) 

13(q) 

Engine  Parameters 

nf 

constant 

- 

(W2E)IC 

constant 

- 

W2EC 

f(N  F/yfT2) 

13  (m) 

cowl.  For  severe  hammer  shocks,  flow  reversal  occurs  within  the  duct,  resulting  in 
reverse  flow  spillage  at  the  cowl.  Knowledge  of  the  pressure  level  behind  the  hammer- 
shock  is  important  when  designing  subsonic  diffusers  since  it  can  be  considerably 
higher  than  freestream  total  pressure.  The  theoretical  approach  to  hammershock 
simulation  used  in  this  model  is  discussed  below. 


Dynamic  Calculations 

The  only  key  area  in  which  this  model  differs  significantly  from  Model  2  is  in  the 
approach  used  to  simulate  NS  travel  upstream  along  the  ramp  during  a  hammershock. 
Detailed  discussion  of  the  upstream  hammershock  model  is  included  in  the  general 
review  of  the  overall  simulation  approach  that  follows. 

The  three-lump  "lumped  parameter"  approach  provides  an  approximate  simula¬ 
tion  of  the  actual  response  of  the  system  to  a  sudden  cutoff  in  engine  demand.  The 
total  pressures  within  each  of  the  three  lumped  volumes  are  calculated  using  the  same 
approach  as  discussed  in  detail  under  the  dynamic  section  for  Model  1.  When  outflow 
from  Vol18  is  suddenly  cut  off,  the  instantaneous  flow  imbalance  across  the  lumped 
volume  (Volpg)  causes  its  pressure  Bp 2  and  total  temperature  Tp2  to  rise.  Instantane¬ 
ous  outflows  at  each  of  the  two  fixed  boundaries  (Stations  14  and  8)  are  computed  by 

•  • 

integrating  the  rates  of  change  of  flows  (Wl4  and  Wg).  Since  these  rates  are  a  function 
of  the  pressure  levels  across  the  fixed  boundaries,  the  rise  in  pressure  in  Vollg 
reduces  instantaneous  outflow  (Wp4)  from  Vol12,  thereby  producing  an  instantaneous 
flow  imbalance  across  Vol12,  and  increasing  its  pressure  PT12  •  Instantaneous  out¬ 
flow  from  the  fixed  boundary  at  Station  8  is  similarly  reduced  by  the  rise  in  PT12  » 
resulting  in  an  increase  in  Volg  pressure  pT6.  The  normal  shock  is  unaffected  by  the 
sudden  cutoff  in  engine  demand  until  the  rise  in  P-pg  forces  it  to  start  traveling  rapidly 
upstream. 


Figure  10(a)  shows  the  NS  located  downstream  of  the  first  ramp.  The  Volg 
upstream  boundary  is  the  NS  which  extends  to  the  critical  flow  streamline.  The  inflow 
(Wx)  based  on  Ay ,  includes  the  subsonic  spillage,  Wgp .  During  the  period  the  NS  is 
downstream  of  the  leading  edge  (i.e.,  X  >  0),  the  critical  flow  areas  Ay  and  P-px^TO 
are  obtained  from  input  data  curves  as  a  function  of  X.  The  variation  in  P-px^TO  as 
the  NS  travels  from  the  cowl  to  the  inlet  leading  edge  is  obtained  by  preselecting  the 
three  known  pressure  ratio  values  for  NS  at  (1)  inlet  leading  edge  (PTX/pTo  =  1), 

(2)  end  of  first  inlet  ramp,  and  (3)  end  of  second  inlet  ramp,  from  the  inlet  input  data 
at  the  given  flight  conditions.  Linear  variation  in  P-px^TO  is  assumeci  between  these 
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I 


(a)  Normal  Shock  Downstream  of  First  Ramp 


(b)  Normal  Shock  Upstream  of  First  Ramp 


Figure  10.  Model  3  -  Super sonlc/Subsonlc  Upstream  Model 
for  Simulating  Hammershock 
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three  values.  The  effective  external  Mach  number  Mx„  upstream  of  the  NS  is  obtained 

?x  (w/TT/PTA)x|.  For  severe  hammer- 


from  the  corresponding  flow  function  F. 


shocks,  flow  reversal  occurs  within  the  subsonic  duct,  usually  before  the  NS  reaches 
the  inlet  leading  edge. 


In  simulating  the  normal  shock  travel  upstream  along  the  ramp,  it  is  assumed 
that  there  is  no  massive  boundary  layer  separation  (i.e.,  AgX  -  Ax).  This  is  accep¬ 
table  as  a  simplifying  assumption,  since  for  a  hammer  shock  transient,  the  primary 
interest  is  in  the  peak  pressures  and  temperatures  experienced  within  the  subsonic 
duct  rather  than  the  accurate  evaluation  of  the  extent  of  the  normal  shock  travel  on  the 
inlet  ramp. 

Figure  10(b)  schematically  shows  the  instantaneous  inflow  and  outflow  through 
Volg  when  the  NS  is  upstream  of  the  inlet  leading  edge  (X  ?  0).  During  simulation,  as 
long  as  X  ^  0,  the  inflow  into  the  NS  is  assumed  to  be  at  /reestream  operating  condi¬ 
tions  (i.e.,  Ax  =  Ag  ;  Mxe  =  Mq  ;  P-px^TO  =  *)•  As  travels  upstream  from 

the  initial  near-the-cowl  position,  the  instantaneous  inflow  into  Volg  can  increase  sig¬ 
nificantly  due  to  flow  reversal  at  Station  8.  The  instantaneous  imbalance  between  this 
inflow  and  the  outflow  (i.e.,  spillage  plus  bleed  flows)  increases  P^,g  ,  until  the  point 
at  which  the  outflow  becomes  instantaneously  greater  than  inflow  due  to  the  continuing 
increase  in  spillage  flow.  At  this  point,  PTg  starts  to  reduce,  thus  relieving  the 
pressure  buildup  within  the  duct. 

The  overall  procedure  for  calculating  normal  shock  (NS)  velocity  X  is  the  same 
as  for  Models  1  and  2.  A  block  diagram  of  the  equations  used  in  formulating  the  ham¬ 
mer  shock  model  is  presented  in  Figure  11.  The  procedure  for  calculating  reverse 
flows  at  the  two  fixed  boundaries  within  the  subsonic  duct  (Stations  8  and  14)  is  the 
same  as  that  used  for  calculation  in  Model  2.  Option  is  included  in  program  logic 
(Figure  11,  Sheet  5)  to  calculate  Prplg  for  hammershocks  severe  enough  to  involve 

reverse  flow  at  the  compressor  face  (W2  <  0).  An  appropriate  value  of  the  engine  stall 
total  temperature  (T^^g)  is  required  as  input  when  simulating  such  a  transient. 


Initial  Conditions  and  Input  Data 

The  calculation  sheet  for  the  initial  conditions  is  presented  in  Figure  12  and  is 
self-explanatory.  Input  data  are  presented  in  Table  in. 
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Figure  11.  Model  3  -  2DE  Three  Lump  Hammershock  Model  Block  Diagram  - 
°  Dynamic  Section:  Sheet  1  of  6 


She* X  I 
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Figure  11  Continued:  Sheet  2  of  6 
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Figure  11  Continued:  Sheet  3  of  6 


Constant 


Figure  11  Continued:  Sheet  4  of  6 


Figure  11  Concluded:  Sheet  6  of  6 
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TABLE  ID.  MODEL  3  -  2DE  THREE  LUMP  HAMMEIiSHOCK  MODEL  INPUT  DATA 


INPUT  PARAMETER 

FUNCTIONAL  RELATION 

FIGURE  NO. 

ISHiliH 

0 

constant 

- 

a 

constant 

- 

H 

constant 

- 

Mo 

constant 

- 

acap 

eonstimt 

- 

ao/acap 

constant 

QR2 

constant 

VW-> 

constant 

- 

^l/p 

constant 

KSP 

constant 

- 

^WRATE/W2E)* 

f(t) 

24(a) 

PTX/PT0 

f(X) 

22(a) 

Ax 

f(X) 

22(c) 

MXI 

f(Fx  *  /R) 

13  (n) 

tostd 

f(H) 

- 

P0STD 

f(H) 

- 

PSY/PTY 

f(Fy) 

13(p) 

Scale 

constant 

- 

Lumped  Volume  Parameters 

(PT6/PTX)IC 

f(Mx) 

13(o) 

A12 

constant 

- 

> 

h-* 

00 

constant 

- 

X8 

constant 

- 

Vo!6 

f(X) 

22(b) 

V0l12 

constant 

V0l18 

constant 

- 

*1 

constant 

- 

*2 

constant 

- 

Z6’  Z12’  Z18 

f(F0>,  f(F12),  f(P18) 

13  (r) 

CV  C12’  C’l8 

f(F6),  f(F12),  f<Fls> 

1 3  (q) 

Engine  Parameters 

^WlC 

constant 

- 

*The  hammershock  input  transient  profile 


SECTION  IV 


SIMULATION  RUNS 


Typical  external  and  internal  transients,  important  with  respect  to  advanced 
fighter  aircraft  air  induction  systems,  vary  considerably  in  the  rates,  total  amplitudes 
and  the  total  time  spans  over  which  they  act.  In  selecting  the  transients  for  the  simu¬ 
lation  runs,  an  attempt  has  been  made  to  cover  a  wide  range  of  input  transients.  This 
section  describes  the  dynamic  response  of  the  full  scale  2DM  and  2DE  air  induction 
systems  to  the  following  set  of  typical  transients. 

2 DM  air  induction  system:  Gust 

Throttle  Burst 
Buzz 

Hammershock 

2DE  air  induction  system:  Afterburner  Blowout 

Angle  of  Attack 
Buzz 

Hammershock 

The  output  data  presented  was  plotted  automatically  at  the  end  of  each  simula¬ 
tion  run  using  IBM  360  on-line  plotting  capability. 

2DM  AIR  INDUCTION  SYSTEM 

Input  data  for  the  program,  developed  for  the  2DM  air  induction  system  model 
described  in  detail  in  Volume  I  of  this  report,  are  presented  in  Figure  13.  Since  the 
full-scale  air  induction  system  is  to  be  simulated  in  the  program,  multipliers  are 
used  throughout  the  program  to  scale  up  the  1/10  scale  model  data  of  Figure  13.  Only 
those  figures  requiring  an  explanation  are  discussed  below. 

Discussion  of  Input  Data 

In  the  triple  oblique  shock  2D  mixed  compression  inlet,  the  local  (ML)  and 
ramp  (M^)  Mach  numbers  are  obtained  from  Figures  13(a)  and  13(b),  respectively.. 

M^  is  presented  as  a  function  of  Mq  and  (  o+  a  ^)  [where  the  first  ramp  angle 
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Freestream  Mach  Number,  Wq 


(a)  ML  =  f  (M0,a) 


Figure  13.  2DM  Air  Induction  System  Input  Data 


Local  Mach  Number,  M 


(b)  MR  =  f  (Ml,  cR2) 


Figure  13  Continued 


E 


2.2 


2. 


Freestream  Mach  Number,  M 

(c)  PTL/PT0  '  f  (M0,  a) 


Figure  13  Continued 


Critical  Mass  Flow  Ratio,  An/A 


<“>VACAp-f<Mo> 


Figure  13  Continued 


Model  Reference  Terminal  Shock  Location,  Xq..™  -  inches 


Total  Bleed  Flow  Ratio 


Ramp  Bleed  Ratio 


Figure  13  Continued 


ao 


14 


<k)  =  f  (X,  aR2) 


Figure  13  Continued 

L 

f 
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(1)  Vol6  =  f  (  aR2,  X) 
Figure  13  Continued 
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Model  Corrected  Engine  Flow,  W \  <5g  '  *  Ibs/sec 


.8  2.2  2.6  3.0  3.4 

Duct  Supersonic  Flow  Function  Parameter,  F^j-  -fR 


(n)  =  f  (FXI  •  /R) 

Figure  13  Continued 
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1st 


Figure  13  Continued 
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0 


.1 


.2  .3  .4 

Flow  Function  Behind  NS,  F^ 

(p>  fe)=  '(Fy> 

Figure  13  Continued 


74 


Lumped  Volume  Coefficient,  Zc  .  Lumped  Volume  Coefficient, 


Flow  Function,  F„ 
o 

(r)  Z6  =  f(F6) 


Figure  13  Concluded 
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<■<  is  fixed  at  10“  ]  and  MR  is  presented  as  a  function  of  the  second  ramp  angle  oR„ 
and  .  The  Mach  number  behind  the  third  oblique  shock  generated  by  the  cowl 
hp  is  calculated  from  Figure  13(b)  by  entering  it  with  the  cowl  ramp  angle  with  respect 
to  local  flow  (<u  R9  4  3°)  and  the  corresponding  Mach  number  M R .  The  total  pressure 
recovery  I’-px'^TO  is  obtained  similarly  in  three  steps  by  entering  Figure  13(c)  to 
obtain  PTL  'P,^  and  entering  Figure  13(d)  to  obtain  pTR/pT1  :in<1  ptxPTR* 

The  critical  mass  flow  ratio  (A^/A^p)  is  presented  in  Figure  13(e)  as  a  func¬ 
tion  of  Mq  and  is  used  to  calculate  the  critical  flow  (WQ)  entering  the  inlet  during 
mixed  mode  operation. 


Figure  13(f)  shows  the  second  ramp  angle  schedule  as  a  function  of  MQ. 

This  schedule  is  followed  to  maintain  inlet  operation  at  optimum  performance  as  MQ 
increases.  The  geometric  throat  area  schedule  is  obtained  from  the  schedule 
and  is  shown  in  Figure  13(g).  The  desired  NS  location  as  a  function  of  MQ  is  presented 
in  Figure  13(h).  The  NS  should  be  as  close  to  the  throat  as  possible  vet  with  enough 
tolerance  to  avoid  being  disgorged  due  to  a  mild  transient. 


Figure  13(i)  gives  the  total  bleed  flow  ratio  W^/Wq  as  a  function  of  M^  for 
(a)  ramp  bleed  inlet  started,  (b)  throat  bleed  with  NS  located  upstream  end  of  slot, 
and  (e)  throat  bleed  with  NS  located  downstream  end  of  slot.  Since  (Id)  and  (c)  specify 
the  total  bleed  flow  when  the  NS  is  located  at  the  two  ends  of  the  slot,  at  any  inter¬ 
mediate  NS  position,  the  upstream  (!¥_„„)  and  downstream  (\VTCV)  bleed  flows  are 

I  id  A  IBY 

computed  as  shown  in  the  sketch  in  Figure  13(i). 


The  effect  of  second  ramp  movement  on  the  boundary  layer  bleed  slot  geometry 
is  shown  in  Figure  13(j).  The  slot  increases  in  length(S)  as  a  increases,  with  the 
ramp  bleed  length  staying  relatively  constant.  The  NS  reference  position,  with  respect 
to  the  slot,  is  superimposed  on  the  figure  in  the  range  (6  degrees  to  12  degrees) 
for  inlet  mixed  mode  operation. 


Figures  13(k)  and  13  (1)  show  the  subsonic  duct  area  and  volume  (Vol^)  distribu¬ 
tion  as  a  function  of  oR2  and  X.  Vol^2  and  Vol^g  are  assumed  to  be  constant  since 
their  boundaries  are  fixed.  Engine  corrected  flow  (W^^,)  as  a  function  of  corrected 
fan  speed  (Np/V"<92)  is  presented  in  Figure  13(m)  for  the  model  inlet. 


Gust  Transient  (Model  1,  Deck  1) 

The  following  example  is  used  to  simulate  the  dynamic  response  of  the  full  scale 
2DM  air  induction  system  to  a  gust  transient. 
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Transient: 


Mach  Number: 


A  typical  gust  (changes  in  freestream 
velocity)  initiated  at  0.2  second  and 
lasting  for  0. 2  second 

3.0 


Altitude:  61, 000  feet 

Angle  of  Attack:  0 


Control  logic  is  provided  in  this  example  for  two  types  of  unstarts:  choked  throat 
unstart  and  the  unstart  caused  by  disgorgement  of  the  normal  shock  (NS)  when  the 
latter  reaches  the  unstable  position  just  upstream  of  the  throat.  The  inlet  throat  area 
(A^)  is  programmed  to  open  if  the  throat  tends  to  choke  and  a  variable  area  bypass  is 
used  to  control  NS  motion.  At  the  initial  condition,  the  bypass  area  is  sized  such  that 
the  NS  is  located  at  the  desired  supercritical  location. 


Discussion  of  Output  Data.  Output  data  for  this  transient  are  presented  in  Figure 
14.  The  imposed  gust  transient  is  represented  by  reduction  in  velocity  starting  at  0. 2 
second  and  lasting  until  0.4  second  [(Figure  14(a)]  .  The  freestream  static  pressure 
(PQ)  is  held  constant  during  the  simulation.  Hence,  this  transient  essentially  involves 
perturbation  in  freestream  total  pressure  (PTQ)  corresponding  to  the  reduction  in  MQ 
as  shown  in  the  same  figure.  The  run  is  continued  to  0. 7  second  to  obtain  complete 
convergence  to  initial  conditions.  As  discussed  before,  two  independent  types  of 
unstarts  can  occur  when  a  gust  transient  is  imposed  on  the  2DM  air  induction  system. 
The  NS  unstart  is  prevented  by  the  bypass  area  (A_,v)  control  and  the  choked  throat 

-D  I 

unstart  is  prevented  by  the  throat  area  (AT)  control.  These  are  discussed  separately 
below. 

The  effect  of  the  gust  transient  on  the  NS  position,  and  the  action  of  the  bypaj  s 
area  Agy  to  prevent  unstarting,  is  discussed  first  using  Figures  14(b)  and  14(c). 
Initially,  the  NS  is  located  [Figure  14(b)]  8.7  inches  downstream  of  the  throat. 
Reduction  in  MQ  at  0.2  second  reduces  the  NS  upstream  Mach  number  (Mx)  along 
with  a  corresponding  reduction  in  total  pressure.  This  reduction  in  and  total 
pressure  causes  the  NS  to  move  upstream  to  a  smaller  area  (A^).  Figure  14(c) 
shows  that  the  bypass  starts  to  open  in  response  to  the  upstream  movement  of  the  NS, 
thus  increasing  the  bypass  outflow  to  the  point  where  its  effect  becomes  strong  enough 
to  change  the  direction  of  NS  movement.  This  change  in  direction  occurs  at  0.275 
second  when  X  =  154.5  inches,  thus  indicating  that  an  unstart  has  been  avoided  since 
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(f)  Pr^/P-TQ  vs  Time 
Figure  14  Concluded 
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the  NS  is  still  downstream  of  the  throat  (X,^  150  inches).  Having  accomplished  its 

task,  the  bypass  then  starts  to  close  and  since  the  NS  is  still  moving  under  the  influence 
of  the  transient,  both  X  and  A^,  tend  to  overshoot  their  initial  locations.  However, 
once  the  gust  is  over  (at  0. 4  second)  the  initial  values  are  reached  very  rapidly.  The 
bypass  control  logic  error  signal  is  a  function  not  only  of  the  position  error  (X^^-X), 
but  also  of  the  rate  at  which  the  deviation  from  X^ p  occurs  and  also  the  rate  at  which 
the  A_,.,  opens  or  closes. 

r>  l 


The  throat  area  (AT4  response  to  avoid  choked  throat  unstarting  is  now  discussed 
using  Figures  14(d)  and  14(e).  Reduction  in  MQ  causes  a  corresponding  reduction  in 
the  total  pressure.  The  result  is  that  the  throat  flow  function  (WT\  /ttx/atptx)> 

initially  corresponding  to  a  throat  Mach  number  of  1.5,  tends  to  increase  towards  the 
choking  value  (0.532).  The  throat  area  control  is  programmed  to  try  and  hold  the 
input  throat  area  schedule  ATgCI.j  [(Figure  13(g)]  .  Since  this  ATgCH  defines  the 
desired  Frp  as  a  function  of  MQ,  any  deviation  from  this  desired  value  introduces  an 
error  signal  (A-pg^H  -A,^)  into  the  system  which  activates  the  throat  area  control.  In 
this  particular  example,  F^  increases  and  A^  [Figure  14(e)]  responds  by  opening 
until  Ft  returns  to  its  initial  value  at  .0.35  second.  Its  task  accomplished,  AT  then 
returns  to  its  own  initial  value. 


Figure  14(f)  shows  the  instantaneous  pressure  recovery,  PT2/PT0,  during  the 
transient.  As  the  NS  moves  upstream  during  the  initial  phase  of  the  transient  (from 
0.2  to  0.3  second),  the  total  pressure  los's  across  the  progressively  weakening  NS 
decreases,  hence,  improving  the  overall  instantaneous  total  pressure  recovery.  As 
mentioned  previously,  however,  the  absolute  value  of  P^  is  reduced  from  the  initial 
condition  throughout  the  transient  due  to  the  overriding  reduction  in  P^.  A  peak 
value  of  P^/P^q  =  0.8  is  reached  at  0.275  second  before  the  recovery  starts  to  drop 
again  during  the  second  phase  of  the  transient  (because  the  NS  changes  direction  under 
bypass  control)  reaching  a  minimum  value  of  0. 767  at  0.375  second  before  returning 
to  its  initial  value. 


Throttle  Burst  Transient  (Model  1,  Deck  1) 

The  following  example  is  used  to  simulate  the  dynamic  response  of  the  full  scale 
2DM  air  induction  system  to  a  throttle  burst  transient. 

Transient:  Throttle  burst  input  as  a  ramp  signature 

with  percent  fan  speed  increase  from  56.5 
to  75.5  in  a  time  period  of  1.5  seconds 
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Mach  Number: 
Altitude: 

Angle  of  Attack: 


2.5 

61,  000  feet 
0 


Discussion  of  Output  Data.  A  ramp  signature  for  the  percent  engine  fan  speed  (N„) 

—  r 

which  starts  at  0.5  second  and  lasts  for  1.5  seconds  is  input  to  simulate  the  throttle 
burst  [Figure  15(a)] .  The  effect  of  this  change  in  speed  can  be  seen  in  the  increase 
in  engine  flow  demand  Wt>E  [Figure  15(b)]  during  the  transient.  The  effect  of  the 
transient  can  best  be  seen  by  observing  its  effect  on  the  NS  position  X  [Figure  15(e)]  . 

Engine  flow  can  be  increased  either  at  a  high  rate  or  a  low  rate.  The  effect  of 
increasing  the  engine  demand  at  a  high  rate  is  that  the  flow  entering  the  lumped  volume 
does  not  instantaneously  match  the  demand  at  the  lumped  volume  exit.  The  resulting 
instantaneous  excess  of  outflow  over  inflow  results  in  a  rapid  pressure  drop  in  the 
lumped  volumes  and,  hence,  initially  rapid  downstream  movement  of  the  NS  if  not 
compensated  for  with  a  rapid  response  bypass  control  system.  On  the  other  hand, 
increasing  engine  demand  at  a  low  rate  (nearly  steady  state)  has  the  effect  of  merely 
relocating  the  NS  progressively  in  a  new  equilibrium  position.  The  effect  that  dominates 
NS  movement  during  a  transient  is,  therefore,  established  by  the  rate  of  the  input 
transient.  In  this  example,  the  input  transient  is  at  a  relatively  low  rate  with  respect 
to  the  dynamic  response  of  the  air  induction  system  and,  therefore,  the  primary  effect 
is  a  nearly  steady  state  downstream  movement  of  NS.  The  downstream  movement  of 
the  NS,  in  turn,  increases  the  NS  strength  and,  hence,  the  total  pressure  loss  across 
it  as  it  travels  downstream.  This  drop  in  PT2  is  presented  in  Figure  15(d). 

At  0.5  second  when  the  transient  is  initiated,  the  NS  is  located  just  downstream 
of  the  throat  [at  161  inches,  Figure  15(e)]  and  starts  to  travel  downstream  under  the 
influence  of  the  transient.  At  the  same  time  ABY  starts  to  close  [Figure  15(f)]  in 
an  effort  to  return  the  NS  to  the  throat.  The  NS,  however,  continues  to  travel  down¬ 
stream  up  to  171  inches  because  the  influence  of  the  input  transient  is  greater  than  the 
influence  of  the  closure  of  A^.  during  the  time  period  of  the  input  transient.  A  more 
rapid  closure  of  AgY  would  reduce  the  maximum  downstream  excursion  of  the  NS. 

After  the  input  transient  is  over,  the  NS  rapidly  returns  to  its  initial  position  which 

is  the  equilibrium  position  corresponding  to  the  final  value  of  A„v. 

J3Y 

The  instantaneous  imbalance  between  inflow  [Wg,  Figure  15(g)]  and  outflow 
l  W14>  Figure  15(h)]  for  Vol^,  and  similarly  between  inflow  (W^)  and  outflow 
[W2,  Figure  15(i)]  for  Vol^,  is  practically  negligible,  indicating  that  the  transient 
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Time  (seconds) 


(c)  TT2  (°R)  vs  Time 


Figure  15.  2DM  Air  Induction  System  Output  Data  - 
(Model  1,  Deck  1);  M  =  2. 5;  H  =  61000  ft; 


Throttle  Burst 

O'  =  0 


Time  (seconds) 
f)  Aqy  (in2)  vs  Time 


Figure  15  Continued 


(8)  WR  (lbs/ sec)  bs  Time 


Figure  15  Concluded 


is  slow  enough  to  be  nearly  steady  state.  The  amplitude  of  the  maximum  instantaneous 
change  for  the  three  flows  W^,  W^,  and  W2  from  their  initial  values  is  also  not  sig¬ 
nificant.  This  is  because  the  relatively  slow  increase  in  engine  flow  demand  is  met 
by  a  corresponding  decrease  in  bypass  area  at  nearly  the  same  rate,  with  minimal 
dynamic  fluctuations  within  the  subsonic  duct. 


Buzz  Transient  (Model  2,  Deck  3) 

The  following  example  is  used  to  simulate  a  buzz  transient  of  the  full  scale  2DM 
air  induction  system. 


Transient: 

Mach  Number: 
Altitude: 

Angle  of  Attack: 


Buzz 

3.0 

70,  000  feet 

0 


Buzz  is  initiated  by  triggering  separation  through  program  logic  with  the  NS 
located  just  upstream  of  the  throat  at  the  initial  condition.  The  bypass  area  sized  for 
this  NS  location  is  then  frozen  for  the  duration  of  the  run  by  inputting  zero  bypass 
gain  in  the  control  logic. 


Discussion  of  Input  Data.  Additional  input  data  curves,  for  the  buzz  model  (Model  2), 
are  presented  in  Figure  16.  These  data  are  required  for  this  model  in  addition  to  the 
basic  input  data  presented  for  the  Standard  model  (Figure  13).. 


The  total  pressure  recovery  upstream  of  the  NS,  P-j^^TO’  is  presented  as  a 
function  of  NS  location  X  in  Figure  16(a),  During  severe  buzz  the  NS  can  have  large 
excursions  extending  nearly  to  the  tip  of  the  inlet  and  resulting  in  a  noticeable  varia¬ 
tion  of  PTX/PTQ.  Values  of  P^/P,^  are  preselected  from  inlet  input  data  when 
the  NS  is  at  (1)  the  inlet  leading  edge,  (2)  end  of  first  ramp,  (3)  the  cowl,  and  (4)  the 
throat  at  the  given  flight  conditions.  Linear  variation  in  PTX/PT0  is  assumed  between 
these  values. 


Volg  is  presented  as  a  function  of  X  in  Figure  16(b).  It  is  based  on  Ay^,  the  latter 
being  the  critical  flow  cross-sectional  area  when  the  NS  is  external  to  the  cowl  (X<XC) 
and  the  duct  area  when  the  NS  is  inside  the  inlet  (X^XC). 

Figure  16(c)  shows  Av  and  Arv  as  a  function  of  X.  It  can  be  seen  that  the  area 

A  kA 

of  separated  flow  (AEX)  increases  rapidly  as  the  NS  moves  upstream  during  the  buzz 
cycle  (as  discussed  previously  in  Section  III,  the  flow  is  fully  attached  during  downstream 
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Total  Pressure  Recovery,  p-rx//F>TO 


l  » 


Normal  Shock  Location,  X,  in.  (M.S.) 


Figure  16.  2DM  Air  Induction  System  Additional 
Input  Data  for  Buzz  Model 
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Figure  16  Continued 
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(c)  Various  Flow  Areas  as  Function  of  X 


(d)  AAET/AT  =  f(X) 


Figure  16  Continued 
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Static  to  Total  Pressure  Ratio 


^  PS6//PT12  “  f  ^F12^ 


1st  Lumped  Volume  Pressure  Ratio,  P^g/Pgg 


(g)  Ft  *  P-rg/Pgg 


Figure  16  Concluded 


movement  of  the  NS).  Figure  16(d)  shows  the  throat  area  blockage  ratio  AAg^,/A^, 
as  a  function  of  X.  Flow  is  assumed  to  be  attached  when  the  NS  is  downstream  of  the 
throat.  After  an  unstart  is  triggered,  blockage  increases  rapidly  until  the  NS  reaches 
the  cowl,  it  continues  to  increase  at  a  reduced  rate  until  the  NS  reaches  the  first 
ramp,  beyond  which  a  constant  value  of  0.93  is  used.  This  maximum  value  of  0.93 
is  based  on  XB-70  simulation  input  data  used  in  Reference  11  which  compares  XB-70 
flight  test  data  with  simulation  predictions  of  Inlet  unstart  and  buzz.  Spillage  low 
Wgp  is  obtained  using  Moeckel's  equation  [Figure  16(ej. 

Figures  16(f)  and  16(g)  are  used  to  calculate  the  throttled  throat  flow  W__(  ~-W0)'. 

E  T  8 

Figure  16(f)  shows  ra*-‘°  as  a  function  of  the  second  lumped  volume  flow 

function  Fp,.  The  pressure  ratio  includes  the  steady  state  friction  loss  for  the  first 
half  of  the  subsonic  duct. 

Discussion  of  Output  Data.  The  output  data  for  the  buzz  transient  are  presented  in 
Figure  17.  The  run  is  terminated  at  0.36  second  after  simulating  two  buzz  cycles. 
There  are  no  controls  imposed  on  NS  movement  during  the  run. 

When  buzz  is  triggered  through  program  logic  at  the  beginning  of  the  simulation 
run,  the  NS  starts  to  move  rapidly  upstream  from  the  initial  location  just  upstream  of 
the  throat  [Figure  17(a)]  .  The  dominating  effect  during  the  unstart  mode  is  the  throat 
throttling  as  the  NS  moves  upstream.  As  the  geometric  throat  area  Ap  is  throttled 
[Figure  17(b)]  the  throat  flow  Wpp  (=Wg)  is  reduced  [Figure  17(d)]  and  induces  an 
instantaneous  two-way  flow  imbalance,  which  increases  Ppg  [Figure  17(g)]  since  out¬ 
flow  (Wg  +  WTBy)  is  instantaneously  less  than  inflow  (Wx)  and  simultaneously  decreases 
PT12  since  outflow  WJ4  [Figure  17(e)]  is  instantaneously  greater  than  inflow  Wg.  As 
the  NS  is  disgorged  from  the  cowl  at  0.018  second,  the  inlet  starts  spilling,  the  throt¬ 
tling  rate  drops  [Figure  17(b)]  and  the  accelerating  NS  increases  the  rate  of  change  of 
Volg.  With  the  Volg  flow  imbalance  now  reversed  [outflow  (Wg  +  W^gy  +  wgp)  instanta¬ 
neously  greater  than  inflow  (Wx)]  and  Volg  expanding  rapidly,  PTg  starts  to  drop.  It 
should  be  noted  that  although  Wg  continues  to  decrease  because  of  throttling,  the  rapid 
increase  in  Wgp  more  than  offsets  this,  thus  increasing  the  instantaneous  net  outflow 
from  Vol  to  a  level  greater  than  Wv. 

O  A 

Instantaneous  static  pressure  ratios  across  the  NS,  PCY  and  P  (which  are 

OA  b  Y 

themselves  functions  of  the  corresponding  instantaneous  total  pressures  across  the 
NS,  PpX  and  Ppg  respectively)  control  the  NS  fluctuations  on  the  ramp.  Figure  17(a) 
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Figure  17  Continued 


Time  (seconds) 

(1)  TT2Yr)  vs  Time 


Figure  17  Concluded 


shows  that  the  NS  fluctuates  once  between  0.035  and  0.06  second  before  attachment 
occurs  and  triggers  the  second  half  of  the  buzz  cycle.  The  duration  of  time  the  flow 
stays  either  attached  or  separated  is  indicated  on  the  figure. 

At  0.06  second  the  key  buzz  parameter  p>j«12^PT0  [Flgure  17  (c)]  reaches  the 
value  of  0.315  when  reattachment  instantaneously  occurs.  With  no  throttling,  outflow 
from  Volg  rushes  into  Vol12  which  is  at  a  much  lower  pressure  pT12  [Figure  17(h)], 
thereby  causing  a  sucking  effect  on  the  NS,  which  reenters  the  inlet  and  surges  to  a 
highly  supercritical  position  (X  =  170  inches  at  t  =  0.09).  It  reverses  itself  at  this 
point  and  travels  upstream  until  it  again  crosses  the  throat,  thus  initiating  another 
buzz  cycle.  The  reason  why  separation  is  not  retriggered  at  PT12^PT0  0,519 
[see  Figures  7(c)  and  17(c)]  is  because  at  that  time  (0.084  second),  the  NS  is  down¬ 
stream  of  the  throat  and  the  flow  is  assumed  to  be  attached. 

The  buzz  frequencies  for  the  full-scale  and  model  inlets  are  7  and  70  cps 
respectively.  The  total  pressures  and  temperatures  at  the  compressor  face  during 
the  buzz  cycle  are  presented  in  Figures  17 (i)  and  17(1),  respectively.  The  high  rate 
of  change  of  these  properties  and  the  large  amplitudes  (11  psi  for  PT2,  250°R  for 
TT2),  associated  with  severe  buzz  conditions,  are  the  primary  reasons  for  causing 
compressor  stall.  The  pressure  and  the  temperature  variations  increase  the  dynamic 
distortion  at  the  compressor  face  to  intolerable  levels  and  drastically  reduce  the 
compressor  stall  margin. 


Hammershock  Transient  (Model  3,  Deck  5) 


The  dynamic  response  of  the  full  scale  2DM  air  induction  system  to  a  hammer- 
shock  transient  has  been  studied  at  the  following  conditions. 


Transient: 

Mach  Number: 
Altitude: 

Angle  of  Attack: 


100  percent  airflow  cutoff  in  0.01  second 
3.0 

70,000  feet 
0 


No  additional  input  data  are  required  to  simulate  a  hammershock  transient.  At 
the  initial  condition,  the  bypass  area  is  sized  such  that  the  NS  is  located  at  the  desired 
supercritical  location  and  then  frozen  at  this  value  for  the  duration  of  the  run.  The 
bypass  control  logic  is  not  included  in  this  model. 
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Discussion  of  Output  Data.  The  output  data  are  presented  in  Figure  18.  The 
input  transient  [Figure  18(a)]  simulates  engine  stall  by  cutting  off  engine 
airflow  (Wgg)  within  ten  milliseconds.  Thus,  the  air  inuuction  system  continues  to 
vent  a  small  amount  of  airflow  [W9,  Figure  18(f)]  through  the  bypass  area  even  after 
engine  demand  is  reduced  to  zero. 

At  0.01  second,  as  the  Vollg  outflow  W2  starts  to  reduce  at  a  rapid  rate,  the 
inflow  W14  [Figure  18(e)]  is  initially  unaffected,  resulting  in  an  instantaneous  unbalance 
across  the  volume,  which  increases  its  pressure  PT2  [Figure  18(i)].  At  this  time  Pj12 
[Figure  18(h)]  and  P^g  [Figure  18(g)]  are  essentially  unaffected.  As  P^2  continues  to 
rise,  starts  to  reduce  (at  0.013  second)  since  its  rate  of  change  is  a  function  of 
the  total  pressures  across  the  fixed  boundary  (Station  14).  This  reduction  in  VolJ2 
outflow  now  causes  an  instantaneous  flow  imbalance  across  Vol12  and  starts  to  increase 
PT12  at  0,015  second*  Similarly,  a  reduction  in  Volg  outflow  Wg  [Figure  18(d)]  at 
0.02  second  starts  increasing  PTg  three  milliseconds  later.  PT2  continues  to  increase 
until  0.04  second  [Figure  1 8 (i )]  and  then  starts  to  reduce  rapidly  when  the  Volg  instan¬ 
taneous  outflow  becomes  greater  than  inflow  due  to  the  continuously  increasing  sub¬ 
sonic  spillage  as  the  NS  moves  upstream  of  the  cowl.  Note  that  the  subsonic  duct 
experiences  pressures  greater  than  freestream  stagnation  pressure  [Figure  18(c)]and 
that,  due  to  the  adiabatic  compression  of  the  air  within  the  duct  to  pressure  levels 
above  the  equilibrium  recovery  level,  instantaneous  temperature  levels  are  significantly 
greater  than  freestream  stagnation  temperature. 

2DE  AIR  INDUCTION  SYSTEM 

Input  data  for  the  program,  developed  from  the  2DE  1/8  scale  wind  tunnel  model 
described  in  Volume  I  of  this  report,  are  presented  in  Figure  19.  Multipliers  are 
used  within  the  program  to  scale  the  inlet  up  to  full  scale.  Only  those  figures  requiring 
an  explanation  are  discussed  below. 

Discussion  of  Input  Data 

Figure  19(a)  is  obtained  from  shock  tables  and  gives  the  Mach  number  behind  the 
first  oblique  shock  (M^)  as  a  function  of  MQ  and  (a  •»  aR^)  where  the  first  ramp  angle 
is  fixed  at  10  degrees.  A  lower  limit  is  defined  for  the  supersonic  model  (dotted 
line)  since  below  this  line  the  flow  is  subsonic  on  the  second  ramp  due  to  detachment  of 
the  second  ramp  oblique  shock  for  the  given  and  the  corresponding  scheduled  second 
ramp  angle  (oR2). 
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Figure  18.  2DM  Air  Induction  System  Output  Data  -  Hammershock  Transient 
(Model  3,  Deck  5)  MQ  *  3:  H  =  70000  ft:  a  =  0 
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Figure  18  Continued 
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Figure  18  Concluded 


Local  Mach  Number 


For  the  given  aR2  vs  MQ _ 

Schedule,  flow  is  subsonic  on 
second  ramp  below  this  line 


1.6  1.8  2.0  2.2  2.4 

Freestream  Mach  Number,  MQ 


(a)  Ml  =  f  (Mq,  a) 


Figure  19.  2DE  Air  Induction  System  Input  Data 


Local  Mach  Number, 

(d)  P-tx/P-tl  =  f  (Ml, 


Figure  19  Continued 


Critical  Mass  Flow  Ratio,  Aa/A 


Freestream  Mach  Number,  MQ 


(®)  Aq/ ACAp  f  (M0,  oi) 
Figure  19  Continued 
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Model  Inlet  Cowl  to  Ramp  Height,  h,  inches 


(g)  hj  f  (aR2> 
Figure  19  Continued 
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Total  Bleed  Flow  Ratio 


(h) 


=  f  (X,  ml) 


Figure  19  Continued 


ftsch  ~  f 


19  Continued 


16 


20 


4  8  12 

Second  Ramp  Angle,  -  degrees 


^  AXI  f  ^aR2’  ^ 
Figure  19  Continued 
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i 


First  Lumped  Volume,  Vol-  -  inch' 


(k)  Volg  =  f  ( a^,  X) 


Figure  19  Concluded 


The  Mach  number  behind  the  second  shock  (M^g)  is  presented  in  Figure  19(b)  as 
a  function  of  Mg  and  The  total  pressure  recovery  P^/Pjo  is  obtained  similarly 
in  two  steps  as  shown  in  Figures  19(c)  and  19(d).  The  critical  mass  flow  ratio  is  pre¬ 
sented  in  Figure  19(e)  as  a  function  of  MQ  and  (<h  10°). 

Within  the  inlet  operating  speeds,  a  variable  angle  second  ramp  is  required  to 
maintain  optimum  performance  as  flight  Mach  number  increases.  Thus,  some  sort 
of  a  control  mechanism  sensitive  to  flight  Mach  number  is  required  to  accomplish  this 
inlet  geometry  variation  because  the  optimum  geometry  angles  follow  a  definite  Mach 
number  schedule.  This  schedule,  which  gives  optimum  values  of  oR2  for  maximum 
total  pressure  recoveiy  (P^/P^)  as  a  function  of  is  shown  in  Figure  19(f). 
Although  for  this  example  a^2  stays  constant  since  is  constant,  this  schedule  for 
qR2  must  be  followed,  say,  in  an  MQ  transient  (i.e.,  gust)  to  avoid  deterioration  in  inlet 
performance.  An  error  signal  ”  °R2-^ is  &enerafced  by  the  ramp  angle  con¬ 

trol  logic  for  any  deviations  from  the  scheduled  value  and  the  control  is  set  up  to  try 
to  hold  that  value  by  reducing  the  error.  The  model  inlet  throat  height  defined  as  the 
perpendicular  distance  between  the  second  ramp  and  the  cowl  lip  (hp  is  shown  in 
Figure  19(g)  as  a  function  of  a^. 

Figure  19(h)  shows  the  changes  in  boundary  layer  bleed  flow  ratio  (Wg/W2)  as  a 
function  of  NS  position  and  Mg.  Ramp  and  throat  slot  bleeds  are  lumped  together  to 
simplify  the  calculations.  Since  the  bleed  flow  rate  is  a  function  of  the  static  pressure 
ratio  at  the  slot,  a  higher  constant  value  is  used  when  the  NS  is  external  to  the  cowl 
(subsonic  bleed)  and  a  lower  constant  value  is  used  when  the  NS  is  downstream  of  the 
slot  (supersonic  bleed).  A  linear  variation  in  bleed  flow  between  the  above  two  constant 
values  is  assumed  when  the  NS  traverses  the  slot. 


The  bypass  area  is  varied  to  control  the  shock  position.  A  schedule  for  the  throat 
flow  function  as  a  function  of  Mg  is  generated  for  this  purpose  [Figure  19(i)] . 

This  flow  function  is  calculated  ^e^tind  the  NS  when  the  shock  is  located  at  or  slightly 
ahead  of  the  throat  at  various  operating  conditions.  This  schedule  essentially  specifies 
the  desired  subsonic  Mach  number  rt  the  throat  at  various  flight  speeds.  An  error 
function  (F^.g^  -  F^.)  is  generated  by  the  bypass  area  control  logic  for  any  deviation 
of  the  throat  flow  function,  Frf  PTg]  from  the  scheduled  value  FTSCH 

when  the  NS  is  at  a  location  other  than  the  throat.  This  signal  tends  to  close  the  bypass 
when  the  NS  tends  to  travel  too  far  downstream  of  the  throat  (highly  supercritical  oper¬ 
ation  resulting  in  unacceptable  total  pressure  losses)  in  an  attempt  to  drive  it  back 
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to  the  throat  and,  for  the  same  reason,  tends  to  open  the  bypass  when  the  NS  tends  to 
travel  too  far  upstream  of  the  throat  (highly  subcritlcal  operation). 


Afterburner  Blowout  Transient  (Model  1,  Deck  2) 

The  dynamic  response  of  the  full-scale  2DE  air  induction  system  has  been  studied 
using  the  following  example. 


Transient: 


Mach  Number: 
Altitude: 

Angle  of  Attack: 


A  typical  afterburner  blowout  transient  (changes  in  engine 
corrected  flow  demand)  is  initiated  at  1.0  second  and  lasts 
till  1.92  seconds. 

1.8 

72,500  feet 
0 


In  this  example  the  no-bypass  no-control  option  of  the  program  is  exercised  to 
simulate  the  2DE  air  induction  system.  This  is  accomplished  by  assuming  zero  bypass 
area  and  neutralizing  the  bypass  and  ramp  angle  (a^2)  controls  by  inputting  zero  gains 
in  the  control  logic.  With  no  bypass,  there  is  no  means  of  controlling  NS  movement 
and  all  the  excess  flow  is  spilled  subsonically  at  the  cowl  lip  with  the  NS  located  at  the 
proper  standoff  distance  from  the  cowl  to  allow  for  the  required  spillage. 


Discussion  of  Output  Data.  The  afterburner  blowout  at  1  second  [Figure  20(a)] 
causes  engine  exit  air  temperature  to  drop,  thus  rapidly  decreasing  pressure  at  engine 
exit  and  hence,  an  engine  flow  increase  of  about  10  percent  before  it  is  reduced  to  the 
original  level  at  1. 9  seconds  by  reduction  in  nozzle  area.  Since  there  is  no  bypass,  the 
engine  demand  W2E  [Figure  20(b)]  is  the  same  as  diffuser  exit  flow  W2  [Figure  20(f)]. 

The  effect  of  the  transient  on  the  NS  is  shown  in  Figure  20(e).  Initially,  the  NS  is 
located  about  3  inches  ahead  of  the  cowl  lip.  The  rapid  increase  in  engine  flow  demand 
at  1.0  second  causes  depletion  of  the  air  in  the  lumped  volumes  (thus  instantaneously 
reducing  their  total  pressures)  and, hence,  a  sucking  effect  on  the  NS  resulting  in  down¬ 
stream  travel.  The  loss  in  duct  pressure  PT2  [Figure  20(d)]  due  to  Hie  initial  sharp 
increase  in  engine  flow  demand  is  followed  by  a  gain  in  pressure  when  the  engine  flow 
demand  starts  to  decrease  at  1. 1  seconds.  During  the  initial  sharp  rise  of  engine  flow 
demand  (10  percent  within  40  milliseconds),  the  non-steady  effect  dominates  since  the 
pressure  PT2  continues  to  drop,  and  the  NS  continues  to  move  downstream  for  a  short 
period  (60  to  80  milliseconds),  even  after  the  Vol1fi  outflow  W0  reverses  at  1.06  seconds. 
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Time  (seconds) 


(f)  W2  (lbs/sec)  vs  Time 


Figure  20  Concluded 


The  subsequent  reduction  in  engine  demand  from  1. 1  to  1.9  seconds  is  com¬ 
paratively  slow.  At  1. 14  seconds  and  106  inches,  the  NS  reverses  direction  and 
rapidly  travels  upstream  until  it  reaches  the  cowl.  The  sudden  change  in  slope  (dx/dt) 
at  the  cowl  (at  1.4  seconds)  occurs  because  upstream  of  this  point  spillage  flow  effects 
are  significant. 

Angle  of  Attack  Transient  (Model  1,  Deck  2) 

The  2DE  air  induction  system  response  to  this  transient  is  simulated  at  the 
following  operating  conditions. 

Transient:  Angle  of  Attack.  Initiated  at  0.05  second  as  a  ramp 

input  at  the  typical  rate  of  24  degrees/  second  and 
lasting  for  0. 25  second. 

Mach  Number:  2. 2 

Altitude:  50, 000  feet 

The  no-bypass  no-controls  option  is  again  used  for  this  example,  along  with 
the  same  input  data  as  used  for  the  previous  example. 

Discussion  of  Output  Data.  The  input  transient  a  changes  through  6  degrees 
during  the  simulation  run  [Figure  21(a)]  .  The  effect  of  this  change  on  the  properties 
PTX  and  Mx  upstream  of  the  NS  is  shown  in  Figure  21(b)  and  21(c)  respectively.  The 
overall  change  in  PT2,  X,  and  TT<>  as  shown  in  Figures  21(d)  to  21(f)  is  quite  small 
Due  to  the  small  overall  change  in  properties  during  the  simulation,  the  fluctuation 
in  their  values  before  initiation  of  and  after  termination  of  the  transient  is  due  to  the 
exaggerated  scale  of  the  plots. 

As  the  angle  of  attack  changes,  Mx  and  P^  reduce  and  the  NS  moves  upstream 
in  response.  This  loss  in  P^x  is  reflected  in  the  P^g  curve  [Figure  21(d)]  .  The 
TT2  level  [Figure  21(f)]  changes  slightly  during  the  transient  but  returns  to  its  original 
value  after  the  transient  is  over  since  the  process  is  adiabatic. 

As  the  NS  travels  upstream,  the  spillage  flow  Wgp  increases  as  shown  in  Figure 
21(g).  The  change  in  Wg  [Figure  21(h)]  is  seen  to  be  directly  proportional  to  changes 
in  PT2  and  T^2  since  it  is  obtained  from  a  constant  value  of  corrected  engine  flow, 

W2EC' 

I 

f 
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Figure  21  Continued 


Buzz  Transient  (Model  2,  Peck  4) 


The  following  example  is  used  to  simulate  the  buzz  transient. 


Transient: 

Mach  Number: 
Altitude: 

Angle  of  Attack: 


Buzz 

2.5 

50,000  feet 
0 


At  the  design  Mach  number  of  2. 5,  the  inlet  operates  slightly  subcritically  with 
all  the  excess  flow  being  spilled  at  the  cowl  lip.  Buzz  is  initiated  by  triggering 
separation  through  program  logic . 


Discussion  of  Input  Data.  Additional  input  data  curves  required  for  this  transient 
are  presented  in  Figure  22.  The  approach  used  in  developing  these  curves  is  similar 
to  that  followed  for  the  2DM  air  induction  system  except  for  theAAgT  curve  [Figure 
22(d)],  where  the  maximum  level  of  throat  area  blcckago  is  reduced  to  70  percent  as 
compared  to  93  percent  at  Mach  3  for  the  2DM  inlet.  Reference  11,  which  com¬ 
pares  XB-70  flight  test  data  with  simulation  predictions  of  inlet  unstart  and  buzz, 
indicates  that  model  and  flight  tests  for  various  inlets,  including  XB-70,  have 
shown  some  variation  in  the  magnitude  of  boundary  layer  separation  as  a  function  of 
the  static  pressure  rise  across  the  NS  and,  hence,  the  Mach  number  near  the  cowl. 
Since  the  latter  decreases  with  decreasing  Mp,  there  is  a  corresponding  reduction  in 
the  separation  level  (and,  hence,  blockage)  for  the  2DE  contract  inlet  operating  at 
Mach  2. 5  as  compared  to  the  2DM  inlet  operating  at  Mach  3.  The  70  percent  level 
is  based  on  XB-70  simulation  input  data  available  in  Reference  11. 


Discussion  of  Output  Data.  The  output  data  for  the  buzz  transient  are  presented 
in  Figure  23.  The  key  parameter  prrl2/pxo  loggers  separation  and  attachment 
is  shown  in  Figure  23(c).  After  the  flow  attaches  at  0.37  second  (when  pT12//pT0  droPs 
below  0.315),  the  NS  [Figure  23(a)]  moves  rapidly  into  the  inlet  to  a  highly  supercritical 
location.  Separation  is  triggered  again  when  the  NS  emerges  from  the  cowl  since  both 
the  conditions  [(1)  that  P-pi^^TO  ^  a^>ove  0*519  and  (2)  NS  be  upstream  of  the  cowl] 
are  satisfied  at  this  point.  The  reason  why  separation  is  not  triggered  at  0. 381  second 
when  pT12/pT0  Initially  rises  to  0.519  (just  after  attachment  has  occurred)  is  because 
the  NS  is  downstream  of  the  cowl  at  that  time. 
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Total  Pressure  Recovery 


(**)  ^tx^to  ^ 

Figure  22.  2DE  Air  Induction  System  Additional 
Input  Data  for  Buzz  Model 
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First  Lumped  Volume,  Vol„,  in. 


<b)  Volg  =  f  (X) 
Figure  22  Continued 
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i 


(c)  PX12/PT0  vs  Time 


Figure  23 .  2DE  Air  Induction  System  Output  Data  -  Buzz  Transient 
(Model  2,  Deck  4)  MQ  =  2.  5;  H  =  50CD0  ft;  a  =  0 
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(e)  W  (lbs/sec)  vs  Time 


A*X 


Figure  23  Continued 


Time  (seconds) 

(1)  TT2  (°R)  vs  Time 


Figure  23  Concluded 


V 


The  NS  upstream  travel  rate  [Figure  23(a)]  for  the  2DE  inlet  is  slower  than  the 
2DM  inlet  and  the  buzz  frequency  is  lower;  2. 5  cps  compared  to  7  cps  for  the  full  scale 
2DM  inlet.  The  total  pressure  P,^  and  total  temperature  T^.^  amplitudes  of  16  psi 
and  1S0°E  during  buzz  are  shown  in  Figure  23(1)  and  23(1)  respectively. 


Hammershock  Transient  (Model  3,  Deck  6) 

The  dynamic  response  of  the  2DE  air  induction  system  to  a  hammershock  transient 
has  been  studied  at  the  following  2. editions. 


Transient: 

Mach  Number: 
Altitude: 

Angle  of  Attack: 


100  percent  airflow  cutoff  in  0.01  second 
2.5 

50,000  feet 
0 


The  basic  input  data  are  die  same  as  that  used  for  the  buzz  transient. 


Discussion  of  Output  Data.  The  output  data  are  presented  in  Figure  24.  Since 
no  bypass  is  assumed,  engine  demand  is  the  same  as  diffuser  supply  (W2  =  W2£)  and 
both  reduce  to  zero  within  10  milliseconds  after  flow  cutoff  at  0. 01  second  [Figures 
24(a)  and  24(f)].  Although  PTJ2  [(Figure  24(h)] ,  as  expected,  lags  behind  PT2  [Figure 
24(1)] ,  by  about  3  milliseconds,  PTg  [Figure  24(g)]  actually  leads  PT12  by  about  2  milli¬ 
seconds.  The  reason  for  this  is  that  in  the  2DE  inlet  input  data  the  boundary  layer 
bleed  curve  [Wg/W2,  Figure  19(h)]  is  input  as  a  function  of  [Figure  24(f)].  The 
resulting  rapid  reduction  in  instantaneous  bleed  outflow  from  Volg  (about  20  pounds/ 
second  at  initiation  of  hammershock)  to  zero  within  10  milliseconds  causes  Vol„  total 

D 

pressure  PTg  to  rise  somewhat,  immediately  after  the  hammershock  is  initiated, 
before  dropping  rapidly  due  to  the  subsequent  massive  increase  in  instantaneous  spillage 
outflow.  The  reason  why  this  small  initial  increase  in  PTg  does  not  occur  for  the  2DM 
air  induction  system  [Figure  18(g)]  is  because  the  2DM  bleed  curve  [Figure  13(1)]  is 
input  as  a  function  of  Wg  which  stays  constant  during  the  hammershock. 

During  the  hammershock  the  su.sonic  duct  experiences  a  maximum  pressure  of 
1.22  times  the  freestream  total  pressure  at  0.033  second  [Figure  24(c)]. 
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Time  (seconds) 


(c)  PT2//pT0  vs  Time 

Figure  24.  2DE  Air  Induction  System  Output  Data  -  Hammers  hock  Transient 
(Model  3,  Deck  6)  =  2.  5;  H  =  50000  ft;  a  =  0 


Figure  24  Continued 


Time  (seconds) 

(1)  TT2  (°R)  vs  Time 


Figure  24  Concluded 


SECTION  V 


HAMMERSHOCK  SIMULATION  COMPARISONS 


Hammershock  simulation  and  test  data  available  on  the  North  American  Rockwell 
(NAR)  XB-70  and  F108  air  induction  systems  are  compared  with  corresponding  simu¬ 
lation  runs,  using  Model  3.  A  full-scale  XB-70  air  induction  system  hammershock 
simulation  at  Mach  3  (Reference  12)  is  first  compared  with  a  corresponding  Model  3 
simulation  run  using  XB-70  data  available  from  NAR  (References  13  and  14).  This  is 
followed  with  a  comparison  of  Mach  3. 1  hammershock  test  data  on  the  NAR  F108  wing 
tunnel  model,  with  Model  3  simulated  output  using  F108  input  data,  available  in  Refer¬ 
ence  15.  The  results  are  discussed  below. 


NAR  XB-70  INLET  SIMULATION  VS  MODEL  3  SIMULATION 


This  comparison,  using  the  full-scale  XB-70  air  induction  system,  is  based  on 
the  following  operating  conditions. 


Transient: 

Mach  Number: 

Freestream  Pressure: 
Freestre*  n  Temperature: 
Angle  of  i  ttack: 


100  percent  airflow  cut-off  in  1  millisecond 
3.0 

.  251  psi 
373. 5°R 
0 


The  output  data  obtained  from  the  Model  3  simulation  run  is  presented  in 
Figure  25.  Figure  25  (a)  shows  the  hammershock  input  transient  which  is  initiated  at 
.01  second  with  100  percent  cutoff  in  engine  airflow  within  1  millisecond. 

Figure  25  (c)  shows  the  comparison  in  air  induction  system  pressure  recovery 
(PT2/PT0)  between  Model  3  and  the  NAR  model.  The  latter  is  based  on  the  "moving 
hammershock"  approach,  with  the  volume  between  the  hammershock  and  the  engine 
face  treated  as  a  single  lump.  The  model  locates  a  hammershock  (HS)  slightly  forward 
of  the  compressor  face  at  initiation  of  the  transient.  The  lumped  volume  between  this 
shock  front  and  the  engine  face  is  the  control  volume,  which  increases  as  the  HS  moves 
upstream.  Upon  initiation  of  the  transient,  the  instantaneous  outflow  from  the  control 
volume  is  drastically  reduced  as  compared  to  the  instantaneous  inflow.  The 
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Figure  25.  Hammershock  Simulation  Comparison  -  NAR  vs  Model  3 
XB-70  Full  Scale  Inlet;  MA  =  3;  PQ  =  .  251  psi;  TA  =  373.  5®  R;  a=  0 
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instantaneous  control  volume  pressure  rises,  and  the  HS  moves  upstream  at  a  velocity 
required  to  satisfy  the  instantaneous  pressures  across  the  HS.  Since  this  is  a  single 
lump  model,  the  control  volume  between  the  HS  and  the  engine  face  is  very  small  at 
initiation  of  the  transient,  and  grows  to  include  the  whole  subsonic  diffuser  when  the 
HS  reaches  the  cowl.  On  the  other  hand,  Model  3  assumes  a  constant  'lumped'  volume 
(Vol1Q)  near  the  compressor  face.  Thus,  initially  when  the  NAR  control  volume  Is  much 
smaller  than  the  constant  Voljg,  the  instantaneous  within  the  NAR  volume  Increases 
much  more  rapidly  than  within  '.he  Model  3  volume,  as  can  be  seen  in  Figure  25(c). 

This  rate,  however,  drops  off  for  the  NAR  model,  since  its  control  volume  continues  to 
increase  during  the  transient,  and  may  include  the  whole  subsonic  duct  near  the  end  of 
the  simulation  run.  Both  curves  level  off  for  about  3  milliseconds  at  .032  second  and 
then  rise  again  to  the  maximum  values  of  1.3  at  .069  second  for  Model  3,  and  to  1.21  at 
.056  second  for  the  NAR  model 

Due  to  unavailability  of  output  data  from  the  NAR  model  (other  than  PT2/PT0), 
comparison  of  parameters  other  than  P-j<2^T0  was  not  P0S8^e* 


NAR  F108  INLET  TEST  VERSUS  MODEL  3  SIMULATION 


This  comparison,  using  the  NAR  F108  inlet  wind  tunnel  model  (Reference  15)  is 
based  on  the  following  example. 


Transient: 

Mach  Number: 

Wind  Tunnel 

Total  Pressure: 
Total  Temperature: 
Angle  of  Attack: 


100  percent  airflow  cutoff  in  2.5  milliseconds 
3.1 


14. 7  psi 
610°R 
0 


Figure  26  shows  the  comparison  between  the  F108  inlet  model  test  and  the  Model 
3  simulation.  The  latter  compares  well  with  the  test  results.  It  should  be  noted,  how¬ 
ever.  that  the  test  is  in  itself  a  simulation,  since  flow  is  completely  cut  off  by  a  plug  and, 
hence.it  does  not  portray  the  response  to  an  actual  engine  stall  transient,  where  initial 
rapid  cutoff  in  compressor  flow  leads  to  flow  reversal  and,  finally,  some  positive  flow 
while  the  engine  windmills. 


The  complete  simulation  run  is  presented  in  Figure  27.  The  transient  is  initi¬ 
ated  at  ,072  second,  with  flow  demand  completely  cut  off  2.5  milliseconds  later 
[  Figure  27(a)] .  The  instantaneous  total  pressures  PTg,  PT12  a°d  PT2  are  shown  in 
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Pressure  Recovery 


0  2  4  6  8  10  12  14  16  18 

Time  ~  Milliseconds 


Figure  26.  Haxnmershock  Test  vs  Simulation  Comparison  - 
NAR  F108  Wind  Tunnel  Model  vs  Model  3;  Mq  =  3. 1;  a  =  0 
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(a)  Input  Transient,  Wot,  (Ibs/sec)  vs  Time 
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(b)  X  (in)  vs  Time 
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Time  (seconds) 

(e)  PT2//PT0  vs  Time 


Figure  27.  Hammershock  Simulation  -  F108  Wind  Tunnel  Model: 

Mq  =  3.1;  a  =  0 
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(e)  W, .  (lbs/ sec)  vs  Time 


Time  (seconds) 

(i)  PT2  (psi)  vs  Time 


Figure  27  Continued 


Figure  27  Concluded 


Figures  27(g),  27(h),  and  27(i).  It  can  be  seen  that  P^,6  starts  to  respond  about  2.2 
milliseconds  after  the  transient  has  been  initiated,  and  since  this  pressure  affects  NS 
movement,  the  latter  also  starts  moving  at  about  the  same  time  [Figure  27(b)] .  The 
oscillations  in  X  and  PTg  near  the  end  of  the  simulation  run  tend  to  damp  out  fairly 
fast.  Final  steady  state  conditions  calculated  by  Model  3  will  indicate  total  pressure 
recovery  behind  the  NS  of  .3  (which  is  the  pressure  recovery  behind  a  normal  shock  at 
Mach  3. 1)  with  the  NS  located  three  to  four  inches  upstream  of  the  inlet  leading  edge. 
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APPENDIX  I 


DERIVATION  OF  THE  LUMPED  VOLUME  TOTAL 
TEMPERATURE  AND  PRESSURE  R'.TE  EQUATIONS 

In  the  'lumped  parameter'  concept,  the  air  induction  system  subsonic  duct  is 
divided  into  a  series  of  lumped  volumes.  The  volume  nearest  to  the  diffuser  throat 
has  a  moving  upstream  boundary  (the  normal  shock)  and  a  fixed  downstream  boundary. 

The  remaining  lumped  volumes  have  fixed  boundaries.  The  instantaneous  total  tem-  __ 
peratures  and  pressures  w:‘hin  each  of  these  volumes  are  obtained  during  the  simula¬ 
tion  run  by  numerically  integrating  the  corresponding  rate  equations.  These  rate 
equations  are  derived  by  applying  the  one-dimensional  unsteady  continuity  and  energy 
equations  in  conjunction  with  the  equation  of  state  to  each  bumped ’  volume.  The  pro¬ 
cedure  adopted  here  in  developing  these  rate  equations  is,  first,  to  derive  from  basic 
principles,  the  one-dimensional  unsteady  energy  equation  for  a  control  volume  with  a 
moving  upstream  boundary  and  a  fixed  downstream  boundary,  and  then  to  apply  the 
'lumped'  volume  consent  to  this  energy  equation  for  developing  the  rate  equations. 

The  final  forms  in  which  these  equations  are  used  in  the  air  induction  system  dynamic 
simulation  computer  program  are  obtained  through  some  algebraic  manipulation. 

Total  Temperature  Rate  Equation 

The  approach  used  in  deriving  the  basic  unsteady  energy  equation  [Equation  (5) 
in  this  Appendix]  for  a  control  volume  with  a  moving  upstream  boundary  (the  normal  shock) 
and  a  fixed  downstream  boundary  follows  Reference  16.  It  is  included  here  in  the  original 
form  except  for  some  changes  in  nomenclature. 

The  time-dependent  energy  equation  states  that  the  rate  of  energy  accumulation 
in  an  adiabatic  control  volume  without  shaft  work  is  equivalent  to  the  net  flux  of  energy 
through  its  boundaries,  lessened  or  increased  by  the  rate  of  work  done  at  the  surfaces 
of  the  volume.  The  figure  below  portrays  the  situation.  The  control  volume  (CV) 
extends  from  a  moving  control  surface  to  a  fixed  end  boundary.  In  this  portrayal, 
translation  of  the  moving  surface  lags  the  fluid  motion.  However,  the  derivation 
applies  to  any  control  surface  speed  relative  to  fluid  speed. 
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Time  =  t  Time  =  t  +  At 


The  energy  contained  by  the  system  of  fluid  nresent  in  the  control  volume  at 

time  t  is  Ej  +  Ejj  .  Since  energy  can  be  adde~  to  an  adiabatic  system  in  the  form  of 
t  t 

work  only,  and  since  energy  must  be  conserved: 

*  ***)  -  ‘*V  =0  (1) 

At  time  t  the  control  volume  contains  the  fluid  encompassed  in  Regions  I  and  n. 
At  t  +  At  the  fluid  Regions  I  and  in  occupy  the  newly  formed  control  volume.  There¬ 
fore,  during  the  same  time  interval  At,  the  energy  change  for  the  control  volume  is: 

“  Ecv/t  *  + 


Combining  Equations  (I)  and  (2)  gives: 


cv- 


-  E 


m 


t  * 


L-t  4-  t.t 


_  dlS 


where 


(3) 


dEcv  = 

The  energy  contained  in  the  control  volume  is  made  up  of  internal  energy  (e), 
kinetic  energy  (u2/2g) ,  and  potential  energy,  the  latter  assumed  to  be  negligible. 


Thus,  from  the  calculus 


where  and  Xq  denote  the  control  volume  entrance  and  exit  stations,  respectively. 

Similarly,  the  energy  flux  across  the  boundaries  is  made  up  of  internal  and 
kinetic  energy,  so  that  the  time  derivative  form  of  Equation  (3)  can  be  written  as: 


±  ji**  uV;.3)fsAdX-fs;A,tU;  .  ^ 

X(t;  *  . 

4-  Ac  U0  C£o  +  u<-  A  3)  “  ^  w/dt  * 

Equation  (5)  is  the  basic  unsteady  energy  equation  for  a  control  volume  with  a  moving 
boundary. 

Now  the  work  performed  on  the  control  volume  is  assumed  to  occur  only  at  the 
upstream  and  downstream  flow  boundaries  and  consists  of  flow  work.  At  the  two 
boundaries  the  fluid  exerts  pressure  Pg  on  the  control  volume  and  enters  the  upstream 

(dX  \ 

ui  ~  dt  /  ft/sec  with  respect  to  it  end  leaves  the  fixed  boundary  at 
uQ  ft/sec.  The  flow  work  performed  on  the  control  volume  is,  therefore, 


bt'iirlcl/'-'fy  (\L) 


Ps. 


A.  (<*■.-£) 


(6) 


The  expression 


is  the  flow  (Wj)  entering  the  control  volume,  thus, 


buu  C%i) 


U/^ 


and,  similarly, 


±W\ 

^  bo unxtoCj  (x^) 


(7) 


(8) 
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Therefore,  the  net  work  performed  on  the  control  volume  Is, 
«lV  *.  U'  _  U 

-rr-  *  -  —  —  wo 


Now  substituting  Equation  (9)  into  Equation  (5)  we  obtain, 


«  &  W.  -  Ul. 

fs.  'fco 


%(.*) 


a  w;  C^si  ♦  (  ^»so  +  Us  )  (10) 

*3  '  2_3  ^ 


where  the  enthalpy  h„  =  e  +  =- 
S  *S 


A  (  (*+^)fs*Jx  *  ^crTrc 

at  Xct) 


—  Ulo  C-p  Tj-o 


where 


h«.  *  H_ 
S  2-3 


«  hr  •  £f>TT 


In  order  to  evaluate  the  integral  of  the  energy  in  the  control  volume  accurately, 
the  axial  distribution  of  fluid  properties  through  the  volume  must  be  known.  This  dis¬ 
tribution  cannot  be  obtained  by  the  present  theory.  To  overcome  this  deficiency,  the 
'lumped'  volume  concept  is  introduced.  In  this  concept  it  is  assumed  that  the  lumped 


instantaneous  flow  properties  correspond  to  the  instantaneous  average  properties  of 
the  entire  control  volume.  Defining  the  control  volume  properties  in  this  manner 
provides  the  equality: 


where  subscript  c  refers  to  the  properties  (i.e.,  energy,  mass)  within  the  'lumped* 
control  volume. 

Thus, 

+  %  )m‘] 


[ChSc  +  ^  mc] 


ci 

Jp  (m<-  Cprr‘  - 

CpTrc)  - 


(13) 


Applying  the  unsteady  continuity  equation  which  states  that  the  mass  accumulation 

fe 

\dt/ 

we  get 


_c  j  within  the  control  volume  is  equal  to  the  net  flux  of  fluid  through  its  boundaries 


-  t4t  _  Ui( 


or 


mt 


j(^t  -uJo)  <4tr. 


(14). 


Substituting  Equation  (14)  in  Equation  (13)  we  get 


,^[crrTJ(^  -u.)4CJ.^c^c) 


(15) 
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instantaneous  flow  properties  correspond  to  the  instantaneous  average  properties  of 
the  entire  control  volume.  Defining  the  control  volume  properties  in  this  manner 
provides  the  equality: 


where  subscript  c  refers  to  the  properties  (i.e.,  energy,  mass)  within  the  'lumped' 
control  volume. 

Thus, 

atl>‘  *  %  )"\]  *  x  i ^  +  r|)m‘  -  ^  m‘] 

=  jf("V 

(13) 


Applying  the  unsteady  continuity  equation  which  states  that  the  mass  accumulation 

within  the  control  volume  is  equal  to  the  net  flux  of  fluid  through  its  boundaries 

we  get 

-  Ult  _  Wc 
dt 

or 

mc  «  j(^t  -  t-'o)  (14). 

Substituting  Equation  (14)  in  Equation  (13)  we  get 


h  +  '  at- [CP ^  dtJ  -  it 
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Then  Equation  (11)  can  be  written  as 


Cp  -  ^cTTeJ.  (16) 


It  is  assumed  that  Cp  is  constant  over  the  range  of  temperatures  under  consideration. 
Now  the  flow  work  term  on  LHS  of  Equation  (16)  is 


RTSc 


)  ■  0 


K.T 


<*T$ 


Sc  & 


dl" 


-  W*)  At  +  £TSc  (Wt  _  W0)# 


(17) 


Substituting  this  in  Equation  (16)  and  performing  the  indicated  derivative  of  the  first 
expression  on  LHS  of  Equation  (16)  we  get 


Cp  [£Js=j\w.  -Ut)  dt  f  TTc  (W;  -Uc)]  -  iTsc 

—  -  Wo)  a:  Cp  (W*.  TTc'  -  (18) 

To  eliminate  the  static  temper  iture  derivative,  it  is  assumed  that 


d  Tsc  c  Tsc  A  TVc. 

d  L'  "fife  d  L" 


(19) 


This  assumption  is  permissible  for  moderate  subsonic  flow  velocities. 

Substituting  Equation  (19)  in  Equation  (18)  and  simplifying,  we  get 

^  Jifc  f  **  dl"  £Cp  —  K.  JSc  ~J  u  CZp  "^Tt  Wo 
dl'  J  Ire.  J 

+  (Wv  -  W0)  (j*.Ts<_  -  ’’tc]  ‘  (2°) 
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Now 


m 


&  *  J\  ^  -  w°  ) 


<it 


fscV< 


Psc  Vc 


(21) 


Substituting  Equation  (21)  into  Equation  (20)  and  rearranging 
cXTyc  a  t 

at 


'T*.  -  UloTVc.  )  +  (^Ui^-Ulc)  [fcTSc  —  C'p7^.cJ 


Let 


tfp[CW;TT,  -  W«TrJ  -  (K.  -  K  )  rTc  (,  -  )] 


c  fii*  r 

cr~tnz.  L 


i  - 


g.  T&c 
CP 
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(22) 
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•  o  -  bT%  ] 


(23) 


and 


f  fic.  ~^~r  c  1 

LPre  ‘  Ts<.  J 


(24) 


and 

C 


C  2 


(25) 


Substituting  Equations  (23),  (24)  and  (25)  into  Equation  (22)  simplifies  it  to 


dTxc 

--  K  Ru,  Tt.  -  wcTr, 

)  ■+  C  Ttc 

*Ur 

evcpTc  L> 

or, 

ATjc 

«■  /  wfc  rT-  -  uc  tTc>  ) 

c  Vc  PT C 

cit 

*  VC  Fr<_ 

(26) 
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Total  Pressure  Rate  Equation 


This  equation  is  derived  by  applying  the  equation  of  state  to  the  fluid  within  the 
control  volume 


PSc  Vc  rSc 


or 


P.Sc  -  &  ^ 


Vc 


(2 


Differentiating  this  with  respect  to  time  t,  we  get 


d  Psc 


_  R. 


r  i /.  d  C^c.  ^>c ) 
i  v<_  - ^ - 


dk 


TSc  dsJc 


*3F 


or, 


dPs<.  ,  £  T  !if  _  mc  Tst  dVc  1  (2 

dt  L  ~3F  vc  dt  dt  J* 


Now  substituting 


the  continuity  equation 


=  Ui^.  -  t*Jc 
dt 


(2 


and  the  assumption 


dP&c. 

dt 


P-Sc-  /  ^  PTfc  \ 
Pt-c  \  ***  / 


(3 


into  Equation  (28)  we  get 


d  Prc  -  Ptc  ^mcTSc  f  .  _i  dTsc  _  _L  d  ~|r 

~dT"  "  v:  I  '  Mc  TSt  dt  "  Vc  ~2Ttr  y 


Now  substituting 


C  Vc. 

‘Sc  c 


*  Psc  Vc 


and  the  assumption 


ciTsc 

dt 


Tst  <dTs, 


TTc  dt 


into  Equation  (31)  and  rearranging  yields, 


d  l^~c 
dt 


-  -  tJ0  \  +  frc_  drTc 

Psc  Vc  '  Trc 


fVc  dVc 
ITT  dT 


(32) 


From  Equation  (24) 


PSc 


Tr, 

Prc  ± 


Substituting  this  in  the  first  term  on  the  RHS  of  Equation  (32)  we  get 


A  Ptc. 

~2 T 
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^UJV—  t«J  0  J  -f 


PTc. 

Ttc 


d  TVc, 

A  t 


PTc  dVc 
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(33) 


Form  of  Total  Temperature  and  Total  Pressure  Rate  Equations  Used  in  Computer 
Program 


Through  algebraic  manipulation  an  alternative  form  of  these  equations,  more 
suitable  for  use  in  the  air  induction  system  dynamic  simulation  computer  program,  is 

T 

ATc 

obtained  as  follows:  Multiplying  Equation  (33)  throughout  by  5 —  and  rearranging,  we 

.  xc 

get, 


2. 

gTr-c 

5tVcPrc 


^  W0  —  Wt  ^  = 


dTi 


Tc 


d  t 


Trc  A  Vc_ 

~VZ  “dT 


„  Tpc.  d 

PT-  dt 


(34) 


dT 


From  this  the  total  temperature  rate 


Tc 


dt 


is, 
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d^Tc 


=  Prc  * 


U-  & 

Cl'  £  Vc 


vj;  rT 


Ul  c  1"t 


Pro'  ^  a,  K 


for  P  in  psi,  T  in  °R,  W  in  lb/ sec,  V  in  cubic  inch,  t  in  seconds  and  R  in  ft  lb/lb  °R. 

Note  that  the  signs  of  the  last  expressions  on  the  RHS  of  the  Equations  (39)  and 

dVc 

(40)  have  been  reversed  since  — is  negative  when  X  is  positive  (NS  movement 
towards  compressor  face). 
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APPENDIX  n 

DERIVATION  OF  THE  EQUATION  FOR  THE  RATE 
OF  CHANGE  OF  FLOW  ACROSS  FIXED 
LUMPED  VOLUME  BOUNDARY 
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APPENDIX  II 


DERIVATION  OF  THE  EQUATION  FOR  THE  RATE 
OF  CHANGE  OF  FLOW  ACROSS  FIXED 
LUMPED  VOLUME  BOUNDARY 


In  the  "lumped  parameter"  concept,  the  air  induction  system  subsonic  duct  is 
divided  into  three  lumped  volumes.  The  volume  nea*  the  throat  (Vol.)  has  a  moving 

D 

upstream  boundary  and  a  fixed  downstream  boundary  and  the  remaining  two  volumes 
(Vol12  and  Vol^g)  have  fixed  boundaries. 

To  calculate  the  rate  of  change  of  flow  rate  (Wg  and  W  at  each  of  the  two  fixed 
boundaries  the  total  subsonic  duct  is  divided  into  two  equal  volumes,  each  volume  con¬ 
taining  a  fixed  boundary  (Stations  8  and  14),  from  the  previous  three  volumes  (Figure  3). 
The  rate  of  change  of  flow  rate  at  these  boundaries  is  then  determined  by  calculating 
the  rate  of  change  of  momentum  within  each  volume  due  to  the  instantaneous  net  imbal¬ 
anced  force  acting  on  it  during  transient  conditions.  For  this  calculation,  both  the 
volumes  are  frozen  at  their  NS-at-throat  value  during  the  transient. 


The  equation  for  the  rate  of  change  of  flow  rate  (Wg)  is  derived  in  this  Appendix. 
Exactly  the  same  procedure  applies  in  deriving  the  rate  of  change  of  flow  rate  (W^). 

The  momentum  control  volume  is  shown  in  Figure  28,  in  which  the  duct  volume  of 
interest  (Vol6  plus  one-half  Vol^ )  is  represented  as  an  equivalent  constant  area 
(Ag)  volume  of  length  t^,  defined  by 

I  _  Volfc  <■  '/i  Vc- 1 12,  (1) 

'  =  Ati 


Applying  the  unsteady  momentum  equation 
with  the  assumption  that  Pgg  and  Pg^2  act 
respectively,' we  get, 

,*• 


to  the  constant  area  momentum  control  volume 
at  the  entrance  and  exit  of  the  control  volume 


(2) 


where  the  steady  state  duct  pressure  loss 


(PTfi  -  PT12)  is  defined  as  APL. 
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For  subsonic  duct  flow,  assuming  that  the  velocity  heads  arc  relatively  small  com¬ 
pared  to  the  static  pressures,  it  is  permissible  to  assume  that  the  entering  and  exit¬ 
ing  momentum  terms  are  equal.  Thus,  Equation  (2)  can  be  written  as. 


A_ 

at 


J  y  dx.  -  [Psfc  -  fsa  -  1 


(3) 


Further,  again  due  to  the  relative  magnitudes  of  the  velocity  heads  relative  to  the 
static  pressures,  the  dynamic  pressures  at  the  entrance  and  exit  stations  are  assumed 
to  be  equal.  Thus,  in  terms  of  PTfi  and  PT12»  (3)  can  be  written  as 


ct(  a 

dtj  3 


dx  s  A, 
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APl 
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(4) 
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hi  order  to  evaluate  the  integral  of  the  LHS  of  Equation  (4)  accurately,  the 
instantaneous  axial  distribution  of  flow  rate  (W)  through  the  volume  must  be  known. 
Since  it  is  not  possible  to  obtain  this  distribution  using  the  current  approach,  to  over¬ 
come  this  problem,  the  'lumped”  volume  concept  is  introduced,  hi  this  concept  it  is 
assumed  that  the  lumped  instantaneous  weight  flow  corresponds  to  the  instantaneous 
average  weight  flow  of  the  entire  momentum  control  volume,  where  the  instantaneous 
average  weight  flow  is 


“« •  i  I  ? 


dX 


or 


[  W  dx 
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lx  We 
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(5) 


Substituting  Equation  (5)  into  Equation  (4),  we  get 

£i  -  -  A&  [«r&  '  PT12_  -  ]  (6) 

^  dt  2.  J 

so  that  the  equation  for  the  rate  of  change  of  flow  rate  (Wg)  across  the  fixed  boundary 
Station  8  is 


^  =  w8  -  j-(  m 

where  the  dynamic  impedance  constant,  1^,  is  defined  as 
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I. 


3Ae 


(8) 


The  equation  for  rate  of  change  of  flow  rate  (W^)  across  the  fixed  boundary  Station 
14,  is  also  obtained  using  the  same  procedure  and  is  given  by, 


A  Uta  _ 
<*tt 

where 
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APPENDIX  m 


MOECKEL’S  MODEL  FOR  SHOCK  STAND-OFF/ SPILLAGE  FLOW  CALCULATIONS 
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APPENDIX  III 


MQECKEL'S  MODEL  FOR  SHOCK  STAND-OFF/SPILLAGE  FLOW  CALCULATIONS 

Moeckel's  continuity  method  (Reference  17)  has  been  used  in  the  dynamic  models 
to  calculate  the  expelled  bow  shock  position  ahead  of  the  inlet  cowl.  Moeckel  derives 
expressions  for  both  plane  and  axially  symmetric  bodies  and  then  applies  the  theory  to 

two-dimensional  and  axisymmetric  inlets.  The  theory  as  applied  to  2D  inlets  is  of 
interest  here  and  is  summarized  below. 


The  above  figure  shows  Moeckel's  representation  of  the  inlet  bow  shock.  The 
subscripts  used  in  the  text  are  defined  below. 

c  Centroid  of  streamtube  passing  through  sonic  line 
m  Freestream  location  of  stagnation  streamline 
S  Sonic  point  of  detached  shock 
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SB  Sonic  point  of  body 
8  Conditions  along  sonic  lino 

0  Freestream  conditions 


The  coordinate  ym  defines  the  location  of  the  stagnation  streamline  that  separates  the 


mass  entering  the  inlet  (y  )  from  that  spilling  subsonically  outside  the  inlet  (yOT) 

m  od 


ym>- 


For  2D  inlets,  the  spillage  is  defined  as  r  =  1  -  (y  /y_  ).  Moeckel's  continuity  method 

m  od 

establishes  a  relationship  between  this  spillage  and  the  NS  stand-off  distance  (L)  using 
the  following  procedure. 


Assumed  Form  of  Detached  Shock  Wave 
This  method  assumes  that 


(a)  the  form  of  the  detached  shock  from  ym  to  the  sonic  point  S  is  a  hyperbola. 
Thus,  for  a  given  MQ,  the  angles  </>g  and  Ag  for  sonic  velocity  behind  the  shock 
are  immediately  known; 


(b)  the  sonic  curve  between  the  detached  shock  and  the  cowl  lip  is  a  straight  line 

/  L 


and  normal  to  the  deflected  streamline  at  S  (i.e.,  n  =  \  ). 

b 


For  the  2D  inlet,  the  following  equation  for  the  shock  location 


%b) 


is  developed 


from  the  hyperbolic  equation  of  the  shock  wave  in  terms  off 


'm 


^sb/’ysb, 


ana  r?. 


(  C  'tan  to  ^  _  C  Ian.  —  "tan >7  m 

Y5B  Ys&  is& 

where 

C  £  ^  tan  r4<e.  -  *  )  (2) 

/  o 

C  is  known  since  the  cotangent  of  the  Mach  angle  0  =  vMn  -  1  and  <j>a  are  known 

u  » 

quantities  for  any  given  freestream  Mach  number.  P  enters  into  the  hyperbolic  shock 
equation  since  typically  detached  shocks  are  asymptotic  to  the  freestream  Mach  lines 
at  large  distances  from  their  foremost  point. 


Application  of  the  Continuity  Equation 

To  determine  the  quantity  yg/yCB,  the  continuity  equation  is  applied  to  the  fluid 
that  passes  through  the  sonic  line  by  using  an  appropriate  average  value  of  the  stagna¬ 
tion  pressure  distribution  immediately  behind  the  hyperbolic  shock.  This  average 
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value  (P  )  is  assumed  to  exist  at  the  streamline  which  represents  the  mass  centroid 
s  c 

of  the  fluid  passing  through  the  sonic  line.  For  plane  flow,  this  streamline  enters 
the  shock  at  yc  =  (yg  -  ym)/2. 


The  continuity  equation  for  2D  flow  between  the  sonic  flow  area  (A  )  and  the 

s 

freestream  area  (AQ)  of  the  stream  tube  can  then  be  written  as 


where  B  is  a  function  of  Ihs  freestream  Mach  number  only. 


Independently,  from  inlet  and  shock  geometry,  A  /A  may  be  written  as 

S  0 

*2  .  Vs  ~  VS6  (4. 

Ao  (Vs-  Y„)“s1 

Now  for  plane  flow,  n  is  assumed  equal  to  Xg  so  that, 

As  =  Vs  -  Vsa  „  6  (5) 

Ao  (Vs  -  Vsfr)  coiAs 

or 


.Vs  .  |  - 

Vs 8  |  —  &  c&s/lg 


(6) 


Substituting  this  value  into  Equation  (1)  gives  the  final  equation  in  terms  of  the  subsonic 


spillage. 

L 

.  (,  - 1-  vc  +bs,AXsV 

(7) 

Vs* 

V.  YS6/V  t  -6ccsXsy 

Equation 

(7)  can  be  written  as 

L 

-  ksp  ( i  -  Ip. ) 

(8) 

Ysft 

V  VSB  > 

where 

C  4-  6  SinXs  ^  f 

-  SccsXs 

(9) 
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In  applying  Equation  (8)  to  the  2DE  anil  20M  inlets  described  in  detail  in  Volume  I 
of  this  report,  it  can  be  written  as 


UJ  Sr* 

^ O 


(10) 


with  the  freestream  Mach  number  as  used  in  Moeckel  now  becoming  the  relevant  local 
Mach  number  upstream  of  the  detached  shock. 

ajpe  used  as  program 
input  for  calculating  either  the  shock  stand-off  distance  or  the  subsonic  spillage. 


Values  of  Kgp  as  a  function  of  local  Mach  number  (M 
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